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ABSTRACT 


In recent years there has been much work in rough surface scatter- 
ing. The greatest part of this has been done in the area of electromagnetics 
because of the greater possible rewards from radar. However, there has 
been much progress in acoustic wave scattering and both areas have bene- 
fited from the similarity of the phenomena. This work is a study of the 
scattering of acoustic waves from randomly rough surfaces and the extent 
of the similarity to the analogous electromagnetics problem. 

A theory is developed using the method of physical optics for the 
scattering of acoustic waves from randomly rough surfaces which are im- 
perfectly reflecting yet homogeneous with respect to material parameters. 
This development is made for a surface which is given by a stationary, 
ergodic process with a Gaussian distribution of surface heights. The ex- 
pression for th e far zone- acoustic power density is derived and the separa- 
tion into scattered and reflected power components is made. For a surface 
sufficiently rough, the reflected component becomes negligibly small and 
there is near- total conversion to scattered power. For this case the average 
differential scattering cross section is computed from the scattered power 
density for the cases of exponential and Gaussian dependence of the auto- 
correlation coefficient of surface heights. 

The theory is compared with the results of omnidirectional measure- 
ments using surfaces constructed to satisfy the assumptions of the physical 
optics method. The statistical parameters of the surfaces were determined 
from sampled data measurements and were used in the calculation of the 
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theoretical values. The comparison of theory with experiment is favorable 
and the theory is concluded to be a valid one. 

The similarity between acoustic and electromagnetic re-radiation 
from surfaces is investigated by comparing analogous expressions for 
the re-radiated fields. The expressions for the electromagnetic field 
are derived using the method of physical optics as are the acoustic 
expressions, and the determination is valid only for surfaces to which 
this method applies. Similarities are found to exist and there are cases 
where an exact acoustic simulation can be made. In other situations 
the vector effects of the electromagnetic re-radiation dominate but some 
aspects of the re-radiation can be modeled acoustically. 



CHAPTER I 


INTRODUCTION 


I. A The Nature of Scattering Problems 

The analysis of the re-radiation of acoustic and electromagnetic 
waves from surfaces is a very broad problem with many different phenomena 
to be studied. The simplest form of re-radiation from a surface is reflec- 
tion from an infinite plane. One of the distinguishing features of this re- 
flection is that the re-radiated energy flows in a single direction which can 
be determined from the well known laws. Similar to reflection is the re- 
radiation from other geometric shapes such as spheres, cones, paraboloids, 
etc. The re-radiation from these surfaces differs, however, from the reflec- 
tion from a plane in that t)ie energy flow is not in a single direction but is 
scattered in many. When the incident wave is re-radiated from a regularly 
rough surface such as a corrugated plane (for instance in the shape of a 
sinusoid) , the same scattering of energy occurs but in general it is expected 
to be more diffuse (on an intuitive basis) than that for a smooth but non- 
planar body. This, however, depends on the size of the structure of the 
regular roughness relative to the wavelength of the incident radiation. 

Lord Rayleigh [ 1896], an early investigator of scattering phenomena, 
developed a criterion for surface smoothness, now known as the "Rayleigh 
criterion, " which provides an approximate measure of the degree to which 
energy is scattered. For a plane wave of wavelength X incident upon the 
corrugated surface shown in figure 1. 1, there is a phase difference 

/\ o — - 4~tP hcos ■&- 
» A 

for the two rays shown. When A ^ TP , there is cancellation in the 
direction of mirror reflection (specular reflection) and energy is consequently 
flowing in other directions. Raleigh made the arbitrary choice of A p = tT/j^ 
to establish his criterion for delineation of rough surfaces from smooth ones 
(A^= O ); consequently, the Rayleigh criterion for a surface to be smooth 
is 

h 

8 cos €r 
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Formulation of the Rayleigh Criterion 
Figure 1. 1 

Surfaces which form a class for which the re-radiation process is 
different from that of the regular surfaces are those generated by random 
processes. The difference lies in the way in which the phase elements 
from all parts of the surface combine to produce the re-radiated power. 

For a randomly generated surface which is sufficiently rough, a linear 
superposition of the powers from the elemental surface areas is possible 
while for the regular surface a phasor addition taking into account the 
phase differences between the elements of the re-radiated field is necessary. 
For a random surface not rough enough to produce diffuse scattering, a 
division can be made between scattered power and reflected power, i. e. , 
power for which the elemental phases are important. This difference in the 
re-radiation processes for these two classes of surfaces necessitates the 
use of different techniques for the description,- both theoretical and 
experimental, of. the radiation phenomena associated with each class. 

For a wave incident upon a general surface the field. re-radiated in 
all directions from the surface must be given as a function of the parameters 
defining the media involved. For randomly rough surfaces this implies a • 
statistical description of the field and also of the media. The problem of 
re-radiation from regularly rough surfaces and geometric shapes is a 
"deterministic " one and a statistical description is not required. 
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I. B The Description of the Completed Research 

The goal of this study was the description, both theoretical and 
experimental, of the re-radiation in all directions (omnidirectional re- 
radiation) of acoustic waves ftom randomly rough surfaces and the deter- 
mination of the applicability of the results of the acoustic study to the 
similar problem of electromagnetics. There have been many investigations 
into this general area of study (see Section III A) but all with the exception 
of a few have been restricted to the case of back scattering. However, to 
make a complete investigation it is necessary to determine the character- 
istics .of the re-radiation from rough surfaces omnidirectionally; by doing 
this, all phenomena are observed and an opportunity is afforded to test 
theoretical work over a wider range of angles. 

There are many types of randomly rough surfaces for which the re- 
radiation characteristics can be investigated. These can roughly be cate- 
gorized in one way by relating the size of the surface irregularity to the 
wavelength of the radiation. The investigation was restricted to surfaces 
with structures large with respect to wavelength and that are described as 
"gently undulating. " For surfaces of this type, the re-radiation can be in- 
vestigated using the method of physical optics of the Kirchhoff method as 
discussed in Section II B. This approach has been used before and success- 
fully. However, there have been questions raised recently about the use of 
the method [Hagfors, 1964] and one of the purposes of this work is to 
determine the validity of this approach. 

The re-radiation of acoustic waves from randomly rough surfaces is 
investigated by determining the re-radiated pressure in the far zone through 
the use of the Kirchhoff method. This is done for a surface which is gener- 
ated by a random process which is taken to be stationary and ergodic with 
a Gaussian distribution of heights, as discussed in Section II C and 
Appendix II. The surface is not restricted to be perfectly reflecting, as is 
usually done, but is allowed to be an imperfect reflector; however, it is 
assumed to be of homogeneous material (see Sections II B and III D). From 
the re-radiated pressure field the power density is calculated and the 
separation between reflected and scattered power mentioned earlier is 
made (see Section III C). The reflected power is shown to decrease 
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to a negligibly low value for a surface sufficiently rough (as (Davies , 1954j 
and other workers have shown). Making a restriction to such surfaces, 
the scattered power is then used to calculate the average differential 
scattering cross section, which is a particularly convenient and meaning- 
ful quantity for the description of scattering (see Section III D). This 
quantity depends upon the form of the autocorrelation function of surface 
heights and is calculated for the cases of exponential and Gaussian 
dependences, i. e. , for autocorrelation functions that take the forms 

r = ex p{“c} r««P{-$r} 

To determine the validity of the theory as developed using the 
method of physical optics, the average differential scattering cross section 
was experimentally determined for surfaces specially constructed to be 
gently undulating. Measurements were made omnidirectionally for'two 
surfaces of different roughnesses and different materials; these results 
are shown graphically in Section III E, and the technique of measurement 
and the equipment used are described in Sections II B and III C. To make 
the experiments especially useful in determining the validity of the theory, 
the statistical parameters of the surfaces were measured and the sample 
probability distributions of surface heights were computed. It was found 
that the processes were close to being Gaussian and stationary with the 
sample correlation coefficients having a functional dependence lying between 
Gaussian and exponential behavior, the actual variation being closer to 
Gaussian (see Section III D). The measured values of the statistical 
parameters were used in the theoretical expression for average differential 
scattering cross section and values of this quantity were computed for 
comparison with experiment. This comparison, which is made in Section' 

III E, is a favorable one and the conclusion is reached that the Kirchhoff 
method is valid as it is used in this work. 

The characteristics of the re-radiation of acoustic waves from 
surfaces are known to be similar in some respects to analogous electro- 
magnetic re-radiation. These similarities have been utilized in acoustic 
simulations of electromagnetic problems such as scattering in turbulent 
media and radar backscattering. .To determine how closely the scalar 
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acoustic waves simulate the re-radiation of the vector electromagnetic 
waves, comparison was made between analogous expressions for the 
two cases. The electromagnetic expressions used in the comparison 
were derived through use of the Kirchhoff method and the comparison 
made is valid only for surfaces which are gently undulating. In making 
the comparison it was found that under certain conditions an exact 
simulation of the re-radiation of electromagnetic waves is possible. 

In other situations , simulations can be made with varying degrees of 
approximation. These results are contained in Chapter IV. 
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CHAPTER II THfe REFLECTION AND SCATTERING 

OF ACOUSTIC WAVES FROM ROUGH SURFACES 

II. A The Re-radiation from a Rough Surface; the Several Methods of 
Calculation 

A calculation of the re-radiation of acoustic (and electromagnetic) 
waves from surfaces which are randomly rough must give the quantities 
which describe the process and relate these to the parameters which define 
it. These quantities are the probability distributions of the amplitude and 
phase of the re-radiated field, the mean and variance of the field and the 
mean of its intensity. This does not provide a complete description in a 
statistical sense but from the viewpoint of application it is sufficient. 

The parameters defining the process are those describing the medium in 
which the source and the observation point are located and those 
describing the randomly rough surface. The medium is taken to be an 
ideal fluid defined by its density p and compressibility K . The set of 
parameters defining the surface is composed of two sub- sets which are 
taken to be independent. The roughness of the surface is given by the 
statistical parameters which describe the random height process. The 
material parameters are the other set which complete the description of 
the surface. These parameters can also have a random variation; however, 
the material is taken to be homogeneous and is defined by specifying 
the Lame constants \ ' and p. For the ideal fluid, the rigidity, p , is 
zero and X. ' = K. 

The acoustic field in an ideal fluid is scalar and is defined 
completely by giving the pressure $ . This is seen from Constant Cl 9541 
who gives the following equations which describe the acoustic field 

pSy = --v0 v-v = -K?^ 

where 'v is the velocity vector which has its time derivative given by 
the gradient of the scalar pressure. From these equations, the wave 
equation for <f> is obtained; this is 
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= pK 


Therefore, the solution to the re-radiation problem is a scalar quantity 
which must satisfy the wave equation; at the bounding surface of the 
fluid the situation is more complicated. A scalar wave propagating in an 
ideal fluid and incident upon the boundary of a solid excites both scalar 
and vector (shear) waves which propagate in the solid, which is taken 
to be infinte. 

The boundary conditions that the solution must satisfy are that 
the normal stresses and the normal velocity are continuous across the 
interface and the tangential stress goes to zero[ Ewing, Jardetsky, and 
Press, 1957 ]]. Since the tangential stresses are zero at the surface, 
"slippage " occurs. 

The acoustic wave incident upon the rough, planar surface, 
z = £ (x, y), (see figure 1. 1) is taken to be a plane wave 



where 1< ( = k -tpcos -+■"£ cos -©■<£. ) is the 

incident wave vector, k is the wave number and P = 
is the radius vector of the coordinate system. Harmonic time variation 
of the form e iwt is assumed. The total field at the observation point 
P is 

i > =4-i + ^f 

where d> is the field re-radiated from the surface and is the quantity 
v r 

to be determined. 

The calculation of the re-radiated field is one which has been done 
many times in many different ways, both for the acoustic case and for the 
more complicated electromagnetic case as well. In general, techniques 



The Defining Geometry 
Figure II. 1 

used for the one case are applicable for the other. These techniques 
can be roughly classed into two areas with, of course, some important 
exceptions. These two classifications are: the Rayleigh method, named 
after Lord Rayleigh, and the physical optics or Kirchhoff method. The 
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several methods of calculations falling into these two broad classifications, 
as well as others, are discussed and compared in a review paper by 
Lysanov[ 19583 and in a book by Beckmann [19633'. The methods which are 
Rayleigh methods and several others are discussed here to show the rela- 
tion between these methods and the method of physical optics which is 
discussed extensively and used in the next section. The methods chosen- 
for discussion and comparison were selected to be illustrative and not 
because of any special relative importance. There are several other methods 
not presented such as those of TwerskyCl951a]and Ament [19 5 63 which are 
important but are sufficiently different in approach so as not to contribute 
materially to this discussion. 

The Rayleigh Method 

The Rayleigh method utilizes a superposition of plane waves to 
express the field re-radiated from an uneven (not necessarily -randomly 
rough) surface.. The technique was introduced by Rayleigh [18963 who was 
concerned with the reflection and scattering of normally incident acoustic 
waves from periodically uneven surfaces. This method was extended to 
angles of incidence other than normal and used by La Casce and Tamarkin 
[19563 for the prediction of results of an experiment in which sinusoids 
were used as theTTeflecting surface for acoustic waves. Meecham [ 1956a3 
used a variational technique to improve the accuracy of the solutions and 
Parker [195 63 has generalized the method for the acoustic case by intro- 
ducing integral transform techniques. Rice [19513, working with electro- 
magnetic waves, extended the method to randomly rough surfaces and 
Lapin Cl 9 643 has done a similar thing for the acoustic case. 

.The rough surface is expressed in a Fourier series 
oo 

v ^’Cx ; y) = ^ exp-^-Xq (Vwx -hny 

m^-oc 

V 

where C5= and _A_ is the wavelength of the fundamental 

component which is considered arbitrarily large. The field re-radiated from 
£(x,y) is 
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k> (v*\ f \r\)z. 


where 



n?*< 

n a > k^ 9 "" 


For each set of integers (m,n), there is a plane wave of amplitude A mn 
traveling in the direction given by the direction cosines 

cos-0 0 *= ow co^-e- 0y = co3e 0i = R e b(m^ ) 

k ' k k 

For a randomly rough surface, P(m,n) and the A mn become statistical 
quantities. At the boundary, £ (x,y), the boundary conditions on the 
total field are applied; from this a set of linear equations in the unknowns 
A^^ is obtained. In general, this set is not solvable and approximations 
must be made. To obtain an approximate solution, it is necessary to 
require that 

i^^y)k|«f [ ayow) | « ( | apsy )| « i 

For the randomly rough surfaces, the approximate solution for $ ^ is used 
to compute the statistical quantities which complete the description of 
the re-radiation. This is done by making use of the statistical properties 
of the P(m,n). 

The Rayleigh- method , as far as the formal technique is concerned, 
is an "exact" one; i.e. it does not require the use of approximate boundary 
conditions (as the physical optics method does). However, the method is 
restricted in its application to surfaces which are gently sloping and 
have deviations small relative to a wavelength. 
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The Integral Equation Method 


The integral equation method makes use of Green's theorem to 
obtain an integral equation for the field on the irradiated surface. The 
kernel of this equation, upon solution, is integrated over the surface to 
obtain the re-radiated field. Lysanov [1956] and MeechamC 1956 D, 
working independently, developed this method for the acoustic case and 
the electromagnetic case when the incident wave is horizontally polarized 
and applied it to surfaces which are periodically uneven in one dimension 
and perfectly reflecting. 

The re-radiated field at the observation point P is written as an 
integral using Green's theorem CStratton, 1941D; this is often called 
the Helmholtz integral and takes the form 


&C p )= -L 

ryC 41 p 


[ f §(P IQ) - 4>(Q) Q)].-n <J S 

J L — 


II. 1 


s 

where Q is the source point on the surface of integration S , G(P j Q) is 
the Green's function for free space, $ (Q) is the total field on the surface, 
Vq is the gradient operator which operates on points at the surface and n 
is the unit »':>rmal of the surface. For three dimensional problems the 
Green's function takes the form 


q(p|Q) = ® 


-Xk|Eo--r| 


1 E- 0 - r l 


II. 2 


R 


where R is the vector from Q to P. The total field at P is the sum of 
the incident and re-radiated fields; this is 

<K P )» & p ) + <Pr(P) 

- * X £f[(SjCP|<?) V<, <(K9) 

-4>(Q)VQ'SCP|Q)]--fidS 


II. 3 
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When the observation point becomes a surface point, eq. (II. 3) becomes 
an integral equation with the unknown being the field on S . When the 
field and its normal derivative on S are known, <£ r (P) is found using 
eq. (II. 1). 

The integral equation is greatly simplified if the surface is one 
for which $ (Q) = 0 , or a free surface. For this case, eq. (II. 3) becomes 





l *s 





ci S 


where r” denotes evaluation on S . The unknown for this equation is 
s 

Vq$ (Q) and solution for this quantity completes the problem, for again 
<£ r (P) can be found using eq. (II. 1). The solution for Vq<£ (Q) can be 
accomplished using a Wiener-Hopf technique [ Morse and Feshbach, 19533 , 
but it is first necessary to approximate the kernel. The approximation, 
as in the case of the Rayleigh method, leads to a restriction in the application 
of the method to moderately smooth surfaces. Meecham shows for the two 
dimensional case that the surface must be one for which 


IfcWl* 1 




Lysanov, quoting work not available, obtains a different condition for the 
same problem; he requires 


1 

a. 








Also, from the above, the method is not applicable to surfaces that are 
not free. 


The Method of Small Perturbations 

The method of small perturbations is a technique by which an 
approximation to the field on the mean surface is found and used in the 
Helmholtz integral which now is an integration over this mean surface. 
Miles Cl 9543 has used this method in studying acoustic re-radiation and 
Bass and Bocharovt 19573 and FeinbergC 19513 have applied it to the 
electromagnetic case. 
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The approximate field on the mean surface is found by making a 
transfer of the boundary conditions to the mean surface from the perturbed 
surface and then performing an iteration of solutions. The transfer of the 
boundary conditions is accomplished by expanding the total field, <j> , in 
a power series about the mean surface z = M = 0 with respect to z 
and then substituting this series into the boundary conditions on z ='^’(Y.y), 
The series for <£ is 

(fi~ ^ •+ J_ * . . 

• 3z^ 2s o 

The solution for the total field is taken to be 


<p=: ^>(°> + <p f 




-f* 1 ’ 


where <b is 0( ) (this is the sum of <h , and its specular reflection from 

J (1 ) 1 1 

the mean surface), ' isO(^ s ), etc. In the iteration of solutions in 
the transferred boundary conditions <f> ^ is found in terms of <j> ^ in 

terms of $ - , etc. 

, The success of this method depends upon the deviations from the 

'mean surface being small; Miles found that it was necessary to require 


I kyx,y)|<< 


9 ¥*■#) | « 
3 y 


« r 


Therefore, this method is restricted to the same class of surfaces as the 
Rayleigh method and the integral equation method (Lysanov's class of 
surfaces is different, however) and there is, perhaps, nothing to recommend 
one of these methods over the other. The essential difference between 
the methods lies in the treatment of the boundary conditions. The Rayleigh 
method uses the exact boundary conditions and an approximation is 
required to obtain a solution. The integral, equation method and the method 
of small perturbations are techniques in which approximate fields on the . 
boundary are. determined and then used in the Helmholtz integral. 
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II. B The Method of Physical Optics and the Consideration of Imperfect 
Reflectivity 

The re-radiation of waves incident upon a rough surface cannot be 
determined by the Rayleigh method or by any of the other methods described 
above when the surface irregularities are not small relative to a wavelength 
(Lysanov's method is an exception). This restriction limited greatly the 
scattering and reflection phenomena that could be described and predicted 
until the development of a method by Brekhovskikh [19513 that allows the- 
treatment of large-scale roughness. This is the method of physical optics, 
often called the Kirchhoff method because of the similarity to Kirchhoff's ' 
postulate in diffraction theory. 

Because of its applicability to large-scale roughness, the 
Kirchhoff method has been used often in the solution of practical problems, 
mostly in the area of electromagnetic scattering with many applications to 
scattering from the moon. An early application was that of Isakovich Cl 952] 
who was the first to apply the method to the scattering of acoustic and 
electromagnetic waves from randomly rough surfaces, Brekhovskikh' s 
treatment being restricted to regularly rough surfaces. Working in the 
area of acoustics, Eckart [1953D used the method to describe scattering 
from the surface of the sea. Davies [19543 , in the first use of the method 
to predict radar return from natural surfaces, obtained a rough description 
of his own data on sea clutter [Davies and McFarlane, 1946 3. Several 
important developments in the method were made by HoffmanC 1955b! whose 
results are used later in this chapter. Use of the method in predicting 
radar return from the earth was made by Hayre and Moore [19613. Applica- 
tions to lunar scattering have been made by Hughes [1961 3 , Hayre Cl 961 3, 
Hagfors [1 9 61 , 19643, Daniels[ 1961 , 19633, FungC 19643, Fung and 
Moore [19643 and BeckmannC 19643 . 

The method of physical optics is based on the assumption that the 
incident field is "locally" reflected at a point on the irregular surface as 
if an infinite plane were tangent to the surface at that point; this is known 
as the "Kirchhoff Approximation. " Clearly, this fundamental assumption 
does not restrict the method to roughness of any particular size relative to 
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a wavelength; it does, however, restrict its use to surfaces which are 
"locally flat" with respect to a wavelength. Brekhovskikh's criterion for 
applicability is 

■4lf E? c co.s © » ~A 

where R c is the radius of curvature of the surface at the point, 0 is the 
local angle of incidence that (— lc^ ) makes with the unit normal n , and X. 
is the wavelength. If the material parameters are allowed to vary, a 
restriction as to their variation must also be made. 

There are several complications of a fundamental nature that arise 
in the use of the Kirchhoff method. One of these is the occurrence of 
multiple scattering, or secondary reflection. The method does not deny 
the presence of multiple scattering but, at present, there is no way to. 
take it into account. All earlier applications have neglected it and this 
is done here. The other methods described above allow multiple scatter- 
ing and tacitly account for it. The second complication is that of the 
over- shadowing of one part of the surface by another and the consequent 
reduction of the irradiated area, particularly at low grazing angles. 
BeckmannC 1964] and Bass and Fuks [19643 have recently investigated 
this. In' the work here, restriction is made to surfaces smooth enough 
to allow overshadowing to be neglected. 

The total field at a point on the irradiated surface is written as 
the sum of the incident field evaluated at this point and the "local" 
reflection of the incident’ field. This is 

#f)= fo+icf = <p + e^ rr (i + S?) n.4 

where X 1 is the "infinite plane" or Fresnel reflection coefficient C Stratton, 
19413 . In general, for a wave incident in an ideal fluid upon a semi-infinite 
solid surface, the reflection coefficient is [Ewing, Jardetsky, and Press, 
19573 

■ -p >*- 4 je l} 


II. 5 
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where c 


_ A/fcVp 

<sir> "0- 


■VI 


K£ 


p[_(V) +yH / ]sin s 


- 1 


e — 



Ki 


a *= co+ - © 


- I 


where the primes refer to the parameters of the solid and G is the local 
. - angle of incidence defined above. The unit normal is 

e _ Zx -~pZy 4-^ 

( + 2 X % Zy~ 

where Z and Z are the partial derivatives of z with respect to x and 
x y 

y, respectively, and z = ^ (x,y). Since ^ (x,y) is a function of position 
on the surface, n is also and, consequently, 0 and X 1 • Therefore, 
the reflection coefficient is a random function of the random variable 
and, through this functional relationship and that of eq. (II. 4), the 
total field on the surface becomes a random variable. This was first 
pointed out by Brekhovskikh and later by Aksenov C 1958U and Kovalev 
and Pozdnyak Cl961^ . 

When the total field on the surface is known, the re-radiated 
field at P is given by the Helmholtz integral; using the three dimensional 
Green's function given earlier, this gives 




K 



II. 6 


where $ (Q) is given by eq. (II. 4). At this point, the analysis is usually 

specialized to a perfectly reflecting surface; i. e. , one for which j P| = 1. 

However, this specialization can be avoided if the dependence of _p upon 

Z and Z is brought out and made use of by expanding XT* in a Taylor series 
x y 
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in Z x and about the mean angle of incidence, a. For most surfaces 
of interest, this is the angle between (-Tt^) and 1c , A similar thing has 
been done for the electromagnetic case by Kovalev and Pozdynyak 1 1 9 6Ll 
and Fung, Moore, and Parkins, C 19653. The expansion is 

~ ■+■ BjxZ* •+ Zy 

"*■ ■+ ^ayZy' + * • * II. 7 


where 


q ivSS (&£) G n =z±f3*J?\ 

32 xhx=ol / \QZy}zx- O) x s(^ 2j f'/Zx = o) etC * 

Zy-O 


z r° 


2y = 0 


The coefficients in the expression are not necessarily sniall. This is 
probably best seen by considering a surface rough in only one dimension. 
In this case the B‘s are approximately the slope and higher derivatives of 
the curve of reflection coefficient against angle of incidence, which is, 
in general, complicated CErgin, 19523 . Using eq. ‘s (II. 4) and (II. 7) in 
eq. (II. 6) and specializing to the far field where k>>-^ gives 
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where =T^k . The surface is now restricted to be an aperture of a 
size small enough to permit the approximation 
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yet large enough to have sufficient variation to define the random process 
of the surface. Substitution of this approximation into eq. (II. 8) gives 
the result 
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where K = , and the variation of R is considered to be important 

only in the exponential. Eq. (II. 9) is the far field expression for the 
re-radiated field which results from the use of the method of physical 
.optics. 

The transition from physical optics to geometric optics, and to 
some techniques of calculation CMuhleman, 1964D, is made by applying 
the method of stationary .phase CWilf, 1962] to eq. (II. 9). ' The integral 
for </> r (P) has a stationary point whenever the exponent 

K’ r = f<*x + hC^y + K-^z — (k^x - tix)x -+■ f k oz _- k. (2 _)z. 

is such that 

= Q 

9 x 9y 


The condition for stationarity at a point is then 
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(^OZ'klZ.) + C^QX — k»x)— c) (koz_ _ k|z.)-f (fkcy ~ {< 1 yj*z o 

3* 

or the two conditions 

9&W = _ ( b>x - U) 

0(02. “ tl2 } 

must be satisfied simultaneously at a point (x,y). Applying these conditions 
to the expressions for n , eq. (II. 6), it is easily shown that the stationary 
points are those points that have a normal with the property that 

n = -J\-% 

These are the points at which the fundamental laws of reflection are 
satisfied relative to the observation point P . For example, in the 
specular direction (k ox = k lx , k Qy = k ly , k Qz + k lz = 0) the conditions 
for a stationary point lead to 

gysy) = ay>sy) = 0 

3y Qy 

•or only those parts -of the surface which have tangents parallel to the 
mean surface (taken to be z = 0) contribute and the reflection coefficient- 
is that for the mean plane taken at the angle a . For the backs catter case 
0? o = -k^ ) , the reflection coefficient is that for normal incidence and only 
those regions which have normals oriented back towards the source make 
a contribution to the return. 

The method of stationary phase is not an exact one and the results 
obtained are only approximate. Wilf shows that the correction term to 
the first approximation is 0(k while the approximation itself is 0(k~^^); 

therefore, a large value of k is needed for a good approximation. To 
apply the method of physical optics it is only required that 

or 

R c cos &» -f- 

k 
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which can be satisfied for small k as long as the radius of curvature is 
large enough [Fung, Moore, and Parkins, 19 65] . Therefore, physical 
optics cannot be supplanted by geometric optics except under certain 
circumstances, and methods which do so must be carefully applied. 


II. C The Average Power and it's Separation into Scatter and Specular 
Components 

The field re-radiated from the rough surface, <f> , is a complex 
random variable which depends upon the random process z=5(x,y). 

The statistical description of the re-radiation is facilitated by a resolution 
of into two component fields: the specular, or coherent, field and 
the scatter, or incoherent , field CTwersky, 1963], CMoore, 19573, 
[Hayre, 1962] . The specular field is defined as 

II. 10 

where the brackets denote the configuration average which is taken over 

an ensemble of surfaces. The term coherent is appropriate since B 

and b are determined experimentally in the same manner as for steady, 

coherent fields. Using this definition of the specular field as a basis, 

the scatter field <fi is defined as 
v s 

4>r - + 4s j ^s- II. 11 

where /<j) ~ 0. The intensities of the two component fields can be 

shown to obey the relation 

= S a + 

where the asterisk denotes the complex conjugate. The average far field 
power density of the re-radiated field is given by 




II. 12a 
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II. 12b 


where Z Q = V K / p which is real for the ideal fluid assumed here. In 

terms of the specular and scatter intensities, this is 


<P>= J-(e%D»-)z 0 


The expression for the far field power density is obtained by using 
eq. (II. 9) in eq. (II. 12b). This calculation is done under the restriction 
that the surface is smooth enough to neglect the second degree terms 
of the surface in the expression for <f> . Using the expression for the - 
differential area 

d $ = + 


and’ the expression for ri given earlier, the power density is, under this 
restriction 
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Using the coefficients defined in Appendix I, this is rewritten as 


P = y 

r a. 



Z K- (r-i r') ft 


^1 ■+ 




5 s' 

■rK 





The average of p is a configuration average, to be done over 
an ensemble of surfaces. The random height process, is assumed' to be 
ergodic and the mean is taken by averaging "along" the process. The 
process is assumed to be a "single" one as opposed to a composite surface 
made up of a sum of sub-processes as discussed by Isakovich Cl 9523 and 
Beckmann [1964] . The problem of averaging is further simplified here by 
making the additional assumptions that the random height process is 
stationary and Gaussian with the parameters being the standard deviation 
of heights, <y , the mean, and the autocorrelation coefficient, r . The 
mean is taken to be zero with no loss in generality. The assumption of 
a Gaussian height distribution is justified by the results of measurements 
made on the experimental surfaces, as described in Chapter III. The 
assumption of stationarity for '^(x l y) is, perhaps, not wholly warranted 
in almost any practical situation but it is necessary to make it here to 
obtain tractable expressions for . 

The average of p is carried out using the work of Hoffman [1955 b] 3 
the pertinent parts of which are detailed in Appendix II. Making use of the 


assumptions of ergodicity and stationarity and the expressions in Appendix II, 
the average of the power density given by eq. (II. 14) is 
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where the change of variables 

^ ^ x v-= y'-y' 

has been made and the aperture has been taken to be square with dimension 
2.1. This expression is still general in so far as no explicit assumption has 
been made regarding the autocorrelation coefficient r , except that it must 
fall off quite rapidly relative, to the dimensions of the aperture in order 
that the random process,"^ (x,y), be defined and the average determined 
with high confidence. It must also, because of the stationarity assumption, 
be a function only of u and v not of x and y; i.e. , the correlation of 
the heights does not depend upon absolute position within the aperture. 

The resolution of the average power density into scatter and specular 
components is completed by determining the specular power density and 
then .using eq. (II. 15) to find the scatter power density. From eq's. (II. 10) 
and (II. 13), the specular power density, Py> S p* is 

^ | ^ II. 16 

Making the same approximations on eq. (II. 9) that were used in obtaining 
eq. (II. 14) and using the results of Appendix II, the average field is 
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Therefore, the specular power density is 
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for a square aperture of dimension 21. Davies C19541 has obtained a similar 
expression. 

The effect of the roughness, as described by the parameters a and 
r, upon the relative magnitudes of the specular and scatter power densities 
is determined using eq. (11.17). The power re-radiated from a smooth 
aperture is 



where the integ ratio n takes place over the halfspace above the aperture 
(in the far field). The specular power re-radiated when the surface is 
rough is given by this same integral with cr/0. The ratio of the latter to 
the former is not greater than 


<2Xp 



TTC T CoS 0- lZ . 




Therefore, as cr increases and the surface becomes rougher, a smaller 
fraction of the incident power is specularly reflected and an increasingly 
larger amount is scattered. When is about 

five or greater, for most practical purposes, there is total conversion to 
scatter power. It is important to note that this conversion depends only 
upon a and not r. 

Clearly, the occurrence of 

<0o> ~ Se xb d£ O 



as the roughness increases denotes the total conversion to scatter- 
power. At this point, the field is given by 
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fir- <<{V> + 
De xd 


as described above. Beckmann [19633 has shown that the amplitude D 
of the scatter field is Rayleigh distributed; i.e. 


§(D) = 


fb e - D 4M^ 

J H 3 - 



D^O 

D^o 


where $ denotes the probability density function, and the phase, d , 
is uniformly distributed from zero to 2 v . The uniformity of the phase 
distribution indicates that the scatter field (now the total field) is given 
completely by its intensity which can be found from eq. (II. 15). 

The -re solution of the total power density into components is 
accomplished in a more straight-forward manner by operating directly 
on eq. (II. 15), although it is necessary to make use of eq. (II. 17). The 
exponential in the integrand of eq. (II. 15) is written using a power series 
as 
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Similar expansions have been made but in a more specialized way by 
Clarke [ 1963a] and[Fung 1964] . Using this expansion in eq. (11.15), 
the first term of the integral becomes 


Zo(£L 

a. '•4tte a 
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The zeroth term in this infinite series of integrals (after interchanging 
summation and integration) is recognized from eq. (II. 17) to be the specular 
power density. The other terms in this summation beginning with n = 1 
and the remaining terms in the integral are scatter power 'density terms. 

The terms in eq. (II. 19) starting with n = 1 comprise what is often called 
the- "quasi-specular " term primarily, it seems, because it does not 
involve the partial derivatives of the surface. 

The expression for specular power density obtained from eq. (II. 18) 
does not contain the autocorrelation coefficient (because it is the n = 0 
term); however, all other terms, the quasi-specular and the partial 
derivative terms, do. Because of this, the integrand of the integral for 
the scatter power density falls off as the distance from a point on the u, 
v surface increases as long as the partial derivatives of the auto- 
correlation coefficient and their products are approximately of the same 
order as the autocorrelation coefficient itself. Assuming this is so, the • 
integration over u and v is therefore a summation of the contributions 
made by the elements of the surface over which there is correlation of 
surface heights. Because r is a function of u and v only, the inte- 
gration over u and v is a function of x and y only through the 
limits. This function, again because of the assumption of stationarity , 
is a constant function of x and y with variation only near the edge - 
of the aperture. This variation constitutes only a small fraction of the 
total because the aperture is assumed to be large relative to the decay 
of the autocorrelation coefficient. Therefore, the integration over x and 
y is a linear superposition of the contributions from all parts of the sur- 
face with the portions from the edges being small providing the aperture 
is large relative to the variations of the surface. This is in contrast to 
the integration for the specular power density from which it is seen that 
the re-radiation from the different regions of the surface must add in 
phase to constitute a return. Also, it is this linear superposition of 
powers from regions of the surface that marks the difference between 
scattering from surfaces that are randomly rough and those which are 
regularly rough. 
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When the conversion to scatter power is complete, the specular 
field being essentially zero, the re-radiation is described completely by 
giving the variation of its intensity or scatter power density as a function 
of the angle of incidence and angles of observation. This variation is 
given in a very meaningful way by introducing the average differential 
scattering cross section C Kerr, 19513 . This quantity is defined as the 
average power scattered per unit solid angle, per unit incident power 
density, per unit area of the mean scattering surface, and is usually denoted 
by a Q . This definition is motivated by the fact that power is scattered on 
a per unit area basis for a sufficiently rough surface, as described above. 
Clearly, the differential scattering cross section cannot describe the 
re-radiation when both specular and scatter components are present; how- 
ever, the definition can be applied to the scatter component only and a 
reflection coefficient introduced for the specular component. 

Assuming the specular component is negligible, the expression 
for o 0 is obtained by applying the definition to eq. (II. 15). In making 
this application it is necessary to make several approximations. These 
are best made after eq. (II. 15) is simplified through a series of partial 
integrations. The first and second derivatives of the surface with respect 
to u and v are written as 
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Substituting eq. 's (II. 20) into eq. (II. 15) and performing the partial 


integrations gives 
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In Appendix III it is shown that because of the assumption regarding the 

decay of' the autocorrelation coefficient and its derivatives made earlier, 

the "integrated out " terms are of the order of the specular term or lower 

and hence negligible. The definition for or is now applied to eq. (II. 21) 

^ + 2 

by dividing by the incident power density, {$ ) Z /2 , multiplying by • 

2 ° 

R , and taking the integrand of the integral over x and y ; this gives 
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where the variation of the integral with x and y is neglected and the 
integration takes place over the aperture. 

This expression for a Q is put in a more tractable form if the- region 
of integration can be made infinite. The difficulty in doing this is that 
the specular component, while small and negligible, is not zero and will 
contribute. The n = 0 term of eq. (11.19) (the specular term) upon 
integration over infinite limits will, in fact, yield a Dirac delta function 
which is a misleading result. The problem can be resolved in two ways. 
One of these is to use the expansion given by eq. (II. 18) and then simply 
neglect the n = 0 term when the region of integration is made infinite. 

The other solution is to approximate the autocorrelation coefficient in 
the exponent by the first several terms of a power series about u = v = 0; 
a similar thing has been done by DanielsC 19613 , HughesC 19623 and 
Fung Cl 9 64 3. The approximation is 
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where A^ and A2 are greater than or equal to zero in order that j r j < 1 

near the origin. Of the two possible courses, the latter is chosen to 
avoid a series solution for <7 . Substitution of eq. (II. 23) into eq. (II. 22) 


gives 


CO 
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This approximation for cr Q is recognized to be the Fourier transform 

(within a constant, depending upon definition) of the integrand with respect 

to K and K . 
x y . 

The solution for cr o is completed when a choice is made for the 
autocorrelation coefficient, the A's evaluated, and the integration of 
e.q. (II. 24) carried out. There are a rather. large number of possible 
autocorrelation coefficients from which to choose that have been success- 
ful in describing stationary rough surfaces [Hayre, 1962] . The choice 
is limited, however, to functions which can be closely fit at the origin, 
u = v = 0 , by a small number of terms and that are of such a nature that- 
the resulting approximation can be integrated analytically. The integration 
of eq. (II. 29) can be done analytically only when Aj of A 2 is zero. The 
functions which give this tractability and at the same time are successful 
in the description of rough surfaces are the Gaussian function 

2 

where A^ = 0 and A 2 = l/L and the exponential function 

where A 2 = 0 and A^ = l/L. The quantity L- is commonly called the 
"correlation distance. " Clearly, any other function which has a linear 
or a parabolic behavior near the origin yields a result under the 
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approximation given by eq. (II. 23) which cannot be distinguished from 

that obtained using the exponential or the Gaussian functions, respectively. 

The calculation of a for the approximation to the Gaussian 

° 2 

function is made by substituting = 0 and A 2 = l/L into eq. (II. 24) 
and then integrating using the tables of Bierens de Haan [jl.939j . The 
result ie 



11.28 

The terms are arranged in the following order; the quasi-specular or 

zero degree term is first, the second and third terms are of the first degree 

arising from the first derivatives of the autocorrelation coefficient, and 

the remaining terms are of the second degree. 

The relative importance of the terms of eq. (II. 28) and the behavior 

of the important exponential factor are determined primarily by the magnitude 

2 

of the quantity 2(cr/L) and the manner in which K , K , and K appear, 

2 x y z 

as well as their variation. The quantity 2(a/L) is the variance of 
surface slopes as they would be measured or. a profile made along any cut 
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of the surface as discussed in Section III D. As the surface becomes 
smoother and the slopes smaller the variance of the slopes decreases. 

This results in the exponential factor decreasing more rapidly as the 
angle of observation is changed from the specular direction, where K - 

X 

Ky = 0 and the exponential is unity. Since all terms except the quasi- 

specular term have either or or both as a factor and are therefore 

zero in the specular direction, the qua si- specular term is the dominant 

one and it takes on the maximum value for cr Q because of the behavior of 

the exponential. As the angle of observation moves away from the specular 

direction, the other terms increase in size and. diminish the qua si- specular 

term in importance. This is so because, roughly, the numerators of 

these terms are increasing and the denominators decreasing as the angle 

moves from the specular direction and decreases from the vertical. 

Finally, as the grazing angle is approached away from the specular 

direction, it is the second degree terms which dominate. This is so to 

2 2 

the extent that K does not decrease to the point that K cr is so small 
z z 

that the approximation upon which the result is based is invalidated. 

The effect of the reflectivity of the surface at any observation 

point is determined by the relative size of the terms of eq. (II. 28). The 
2 

constant a^ of the quasi- specular term involves only the value of the 
reflection coefficient evaluated at the angle of incidence measured with 
respect to the mean normal It. The constants of all other terms are 
determined in part by the rate of variation of the reflection coefficient 
with angle of incidence. This means that the power scattered in the 
specular direction is diminished by a factor determined only by the 
reflection coefficient, and the magnitude of the reflection coefficient can 
be measured in a quite straight-forward manner. As the observation point 
moves away from the specular direction, the other terms increase in 
importance and the reflection coefficient, of itself, is no longer the 
determining factor. 

It is interesting as well as important to compare the result for 
the Gaussian correlation function obtained through partial integration to 
that result obtained through a direct integration of eq. (II. 20) after the 
approximation of eq. (11.23) is made. Making this approximation, and 
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‘2 

then substituting = 0, k^ = l/L , and eq. (II. 25) into the approximated 
integral and performing the integration gives the result (a- straight-forward 
but tedious process) 
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Examination of eq. (II. 29) shows that as the phase modulation, K a , 

z 

increases and becomes large with respect to unity this result tends toward 

eq. (II. 28). This is as expected because for the "integrated out" term 

2 2 

to be negligible it was necessary to require exp(-K <y ) to become small. 

z 

It is recalled that this is also a necessary condition for the calculation of 
the scattering cross section. 

The calculation of a Q for the approximation to the exponential 
correlation function proceeds similarly to that for the Gaussian function. 
Substitution of A 2 = 0 and = l/L into eq. (II. 24) and the integrating 
using the tables of Bierens de Haan [19393 gives 
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where the arrangement of the terms is like that of eq. (II. 28). 
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The behavior here is similar to that obtained using the Gaussian 
function in that the quasi- specular term dominates in the specular direction 
and then diminishes in importance at side and backscatter with the 
attendant effects due to the reflectivity of the surface. The quantity 
which determines the rapidity with which cr 0 falls off as the angle of 
observation changes from the specular direction is not the variance of 
slopes but is instead. 



This, perhaps, can be considered to be indicative of the roughness. in 
two ways: one of these is through the "slope" <y/L , and the other is 
through the phase modulation K z <y . 

A very tangible difference between the two results, eq. 's (II. 28) 
and (II. 30), is in the behavior in the specular direction. The Gaussian 
result has no variation other than that due to the change in reflection 
coefficient with angle of incidence and for a perfectly reflecting surface 
would be a constant. This contrasts sharply with the result obtained 
using the approximation to the exponential function which has the variation 
of the reflection coefficient superposed on the a Q for a perfectly reflecting 
surfaces which increases with decreasing grazing angle. 
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CHAPTER III 


THE MEASUREMENT OF THE SCATTERING 
CROSS SECTION AND COMPARISON WITH THEORY 

HI. A Historical Development of Measurement Techniques 

The re-radiation of waves incident upon statistically, rough surfaces 
is a problem that can be treated analytically in only a few cases. This 
is true principally for two reasons: Firstly, the present methods of 
analysis are applicable generally only when the wavelength is very large 
or very small relative to the size of the irregularities, as discussed 
earlier. Secondly, most of the naturally occurring surfaces are so complex 
that for many applications they defy analytical description, as discussed 
in Section III D. Because of these theoretical difficulties and the urgent 
need for answers to pressing practical problems, many past and present 
experimental programs have been undertaken. The vast majority of these 
programs, past and present, are to determine the characteristics of 
backscattered radar signals from a variety of rough surfaces. Despite 
this preponderance of radar backs cattering experiments, there are few 
if any forms of scattering that have not received some strong experimental 
effort. 

A period of much progress, both in terms of results and measurements 
techniques developed, began during World War II with a series of experir 
ments performed to determine the radar back scattering characteristics of 
the sea. The initial experiments of this series were done by the Radiation 
Laboratory of the Massachusetts Institute of Technology [Kerr, 1951] CClapp, 
19463 [Cowan, 1946J and the Telecommunications Establishment of the 
United Kingdom [MacLusky and Davies , 19453 C Davies and McFarlane, 

19463 . Later measurements over the sea were made by Wiltse, Shlesinger, 
and Johnson Cl 9573 , Macdonald C1956] and Campbell [1958, 1959] . Grant 
and Yaplee [19573 have investigated backscattering from the surface of a 
river. Experimental programs now in progress are being conducted jointly 
by the Naval Research Laboratory and the Applied Physics Laboratory of 
the Johns Hopkins University. 
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Following closely after the initial sea clutter experiments were 
the beginnings of large scale experiments to investigate the radar back- 
scattering from the generally more complicated terrestrial, lunar, and 
planetary surfaces. The terrestrial measurement programs have been 
largely devoted to obtaining data from a wide variety of naturally occuring 
surfaces. Programs were established early by the Radiation Laboratory 
(Clapp, 1946] , the Sandia Corporation (Edison, Moore, and Warner, I960] , 
the Goodyear Aircraft Corporation (Reitz, 1958j _ , the Ohio State University 
[Cosgriff, Peake, and Taylor, 195^ , the General Precision Laboratory 

(Campbell, 1958, 1959] and the Philco Corporation (no reference available). 
Later programs were established by the Ryan Aircraft Company, the U. S. 
Army, and more recently by the University of Kansas. Radar investigations 
of the moon have been made with high frequency since the end of World 
War II. The earliest lunar explorations were made by an organization in 
Hungary (Bay, 1946] and the U. S. Army (Moffensen, 1946] [Dewitt and 
Stodola, 1949] . Since these, investigations have been made by groups 
in Australia [Kerr, Shain, and Higgins, 1949) (kerr and Shain, 195l] , 
the United Kingdom [Browne et al. , 195 6] [Evans, 1957] [Evans et al. , 
1959] [Hey and Hughes , 1959] [Hughes, I960] , and in this country at 
the Lincoln Laboratory of the Massachusetts Institute of Technology 
[Pettengill, I960] [Evans and Pettengill, 1963] , the U. S. Naval Research 
Laboratory [Trexler, 1958] [Yaplee et al. , 1958, 1959] , University of 
Texas in conjunction with the Royal Radar Establishment [Straiton and - 
Tolbert, I960] , and the Jet Propulsion Laboratory. Radar contact with 
Venus has been made by the Lincoln Laboratory [Price et al. , 1959] and 
the beginnings of an exploration program carried out by the Jet Propulsion 
Laboratory (Victor et al. , 196l] . 

Although most radar experiments have been restricted to the back- 
scattering mode for reasons of practical application or because of necessity 
(lunar and planetary experiments), a number of investigations of the forward 
re-radiation mode (reflection and scatter) have been directed toward the 
solution of radio communication problems. A much smaller number have 
attacked the general problem of obtaining a complete description of the 
re-radiation from rough surfaces by making omnidirectional measurements 
of reflection and scatter. The problem of radio communication over a 
rough earth is an old one that has received much attention, both theoretical 
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and experimental, since before World War. II. A discussion of the experi- 
mental techniques and results is contained in the book by Beckmann and 
Spizzichino[l963] . The general problem of radar scatter and reflection 
in all directions is one that has been dealt with experimentally only 
recently. Early experiments were done in the United Kingdom by Manton 
Cl 958] who made measurements of radar signals re-radiated from rippled 
water in a laboratory tank. Experiments carried out over a variety of 
surfaces have been made and are continuing to be made in a program at 
the Ohio State University C Taylor, 1964]. 

Paralleling somewhat the progress made in radar measurement 
programs has been the progress made in various programs undertaken, 
mostly by the U. S. Navy, to investigate by acoustical means the natural 
surfaces which are the boundaries of the ocean. These programs have not 
been nearly as numerous nor as ambitious as the analogous radar programs, 
although often the needs of the two sets of programs have been similar as 
are the measurement techniques themselves. Much of the work done has 
been directed toward determining the characteristics of the backscatter, 
or reverberation, from the ocean bottom at sonar frequencies (100 to 
40, 000 cps). Early work which took the form of backscattering from 
harbor bottoms was done by the Naval Research LaboratoryCUrick , 1953]. 
Later work has been done by Jones et al.C 1963] , Mackenzie Cl 9 61] and 
McKinney and Anderson Cl 9 63] . Of these, the latter two have obtained a 
catalog of returns from a variety of bottoms. Nolle et al.C 1963] have 
performed experiments under laboratory condition with sands of known density 
as the re-radiating bottoms. A smaller number of experiments have been 
made on the surface of the ocean. Urick and Hoover Cl 95 6] obtained 
backscattering data as a function of wind velocity and Mellen [19 64] 
and Lieberman [1 963] have performed measurements on doppler-shifted 
backscattered signals to investigate the spectral density of the wave. motion. 
Additional data have come as a by-product of the acoustical simulation of 
radar return. The simulation program at the University of New Mexico 
[Edison, 1961] has given much information on the statistics of the signals 
backscattered from rough surfaces and also on the return from a variety 
of surfaces by way of determining materials suitable for modeling. 
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III. B The Technique of Measurement 

An experimental investigation of the re-radiation of waves, acoustic 
or electromagnetic, incident upon rough surfaces requires,in general, the 
description of both the reflection and scatter components of the re-radiated 
field. This is done through a determination of the reflection coefficient 
and the average differential scattering cross section for all angles of 
incidence and observation as a function of the frequency of the' incident 
radiation. For the case of electromagnetic waves, the polarization of 
the incident and re-radiated fields must also be considered, as discussed 
subsequently. The case considered here is that of re-radiation from a 
surface sufficiently rough at the frequency used that the amount of energy 
reflected is negligible. Therefore, it is sufficient to determine only the 
scattering cross section. This is done as a function of the angles of 
incidence and observation (see figure III. 1) through performance of experi- 
ments based upon the equation derived from the definition of the average 
differential scattering cross section. This equation, which is commonly 
called the radar equation, or, depending upon application, the sonar 
equation, is for the general case of bi-static scattering considered here 

<**#>, 4,)>= f f J A m , 

y J fV*" 4-rr 

S 

where: P($ o , $ $ 2 ^ * s t * ie avera 9 e °f the received power obtained 

over an ensemble of surfaces. 

P T is the transmitted power. 

G,p is the gain function of the transmitting antenna (or 

transducer) which is assumed to be circularly symmetric 
and a function of ^ only. 

G r 0$^) is the gain function of the receiving antenna, also 
assumed to be circularly symmetric. 

\ is the wavelength of the radiation, 

S is the mean surface upon which the radiation is incident . 
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. The integrand of eq. (III. 1) is written in the following order: 
the first term is the incident power density, the second term is the 
scattering cross section, and the third term is the receiving antenna 
aperture divided by the square of the distance, which is the solid angle 
subtended by the receiving aperture. The incident power density is that 
of the spherical wave front at the surface as modified by the gain function 
of the transmitting antenna. It is recalled that the scattering cross section 
is defined on the basis of a plane incident wave front, primarily for 
mathematical simplicity. As long as the region over which there is 
correlation of surface heights is small enough so that the spherical wave 
front is approximately plane the difference is negligible. 

It is perhaps a misnomer to call eq. (III. 1) the radar equation 

% 

because the vector nature of the electromagnetic scattering process is 
not evident in this equation. However, this phraseology has been used 
often in the past, and it has only been fairly recently that vector scatter- 
ing has received much attention CThe Ohio State University, 19 63]. The 
vector nature of the scattering is brought out by introducing separate scatter- 
ing cross sections for the description of the fields scattered in and normal 
to the plane of incidence for each of the two orthogonal polarization com- 
ponents. Clearly, there is no difficulty encountered in calling eq. (III. i) the 
/ 

sonar equation because of the scalar nature of the acoustic scattering process 

Bi- static measurements of the "differential scattering cross section 
can be made for both acoustic and electromagnetic waves through the use 
of a transmitting or receiving- antenna (can be both) that has a beam width 
(between half power points) only several degrees wide. This enables the 
approximation to be made that the scattering cross section is constant 
over the area defined by the intersection of the cones of the antenna 
patterns and the mean surface S , which can be called the "effective 
illuminated area. " For the bi -static measurement considered here there 
is no time variation of the amplitude of the incident radiation over the 
effective illuminated area and the steady state is assumed to exist. This 
is shown in eq. (III. 1) by treating Rj as a constant. This existence of 
the steady state in pulsed radar or sonar work is referred to as "beam 
width limitation. " This is to be contrasted with "pulse width limitation " 
in which the area defined on the ground by the antenna gain functions is 



43 



Geometry of Bi- Static Scattering 
Figure III. 1 ' 
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not completely illuminated at any instant of time due to the shortness 
of the physical extent of the radar of sonar pulse jMoore and Williams, 
1957] . 

Making the restriction of high directivity, eq. (III. 1) becomes 

Qr^ r ) 

S 

hi. 2 


<£C<fc,<A,#J>= Oo(^^ fa) 



It is convenient to normalize the average received power with respect to 
the power received upon direct transmission of the power Prj-.. This 
power is 


T> _ P T G^O) q r (o) -X 

id ~ — , 

(4-ir r n ) 


m. 3 


where r is the distance of separation between the transmitting and 
n 

receiving antennas and G^, (0) and (0) are the gains of these antennas, 
respectively. Dividing eq. (III. 2) by the normalizing power, P D , gives 

Pd 



III. 4 

where g^C&p) and g r (^ r ) are the normalized gain functions of the antennas 
and have a maximum value of unity. The integration accounts for the so- 
called "aperture effect" which is the processing of the scattered signals 
by the gain functions of the antennas. 
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The determination of o- Q for a set of angles ($ Q , , $2 ) is 

then a measurement of the average, normalized power and an evaluation 
of the aperture effect at this set of angles. The aperture effect evaluation, 
is a straight-forward calculation which is given in Appendix IV. The 
measurement of the average, normalized power is not so straight-forward 
and it is necessary to carefully design the experiment so that an accurate 
measurement results. It is critically important that the effective area 
illuminated be large enough to permit the surface to have a sufficient 
number of variations to define the random process. ClarkeCl963a^] has 
shown that the dimensions of this area should be at least ten times the 
correlation distance of the random process of the surface. Therefore, 
it is necessary as the beam width of either the transmitting or receiving 
antenna is decreased to improve the approximation of eq. (III. 2) to 
move the antenna back from the mean surface thereby increasing the 
area of the intersection between the mean surface and the cones defined 
by the half power points of the beams. 


III. C The Description of the Experimental Facility and the Conduct 

of the Experiment 

The aim of the experimental investigation is to obtain data with 
which the theoretical results of Chapter II can be compared. The experiments 
consist of the measurement of the average differential scattering cross 
section of surfaces with known statistical properties and surface parameters. 
A complete set of measurements to achieve this goal would be the deter- 
mination of cr o for all values of the set of angles (<j> Q , <j) ± as a 

function of wavelength for a wide range of statistical parameters and 
surface properties. This, of course, is very ambitious and a compromise 
was reached which was designed to bring out features of omnidirectional 
scattering considered to be most important and at the same time provide 
experimental results that could be considered predictable' from' the assump- 
tions of the theory of Chapter II. It was decided to set aside the important 
question of frequency dependence and measure at a single frequency the 
of two surfaces that are smooth enough at this frequency to be called 
locally flat. The surfaces were chosen to be highly reflecting, different 
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in their statistical properties, and, at the frequency used, sufficiently 
rough so that reflection is negligible. The description of the surfaces 
and the measurement of their properties are given in Section III D. The 
series of angles at which a Q is measured is a compromise between 
complete description and undue work. 

The measurement of <j 0 , as discussed earlier, requires the 
measurement of the average of the received power over an ensemble of 
surfaces and the measurement of the antenna gain functions of the 
transducers for the evaluation of the aperture effect. These measurements 
were performed at the acoustic facility of the Remote Sensing Laboratory 
of the Center for Research, Inc. This facility consists of a system of 
two water tanks (one serving as a reservoir) with the associated pumping 
system; a mechanical system for positioning transducers and providing 
motion of the rough surface targets; and an electronic system for the 
generation, reception and processing of ultrasonic signals. The 
arrangement of the equipment about the tanks is shown in figure III. 3 
and the tanks themselves in figure III. 2. 

The mechanical' system shown schematically with laboratory 
parameters in figure III. 4 and photographically in figures III. 5 and III. 6 
is really two independent systems. One of these is the equipment used 
to align and position the transducers in the tank. The operation of this 
equipment is crucially important as the receiving and transmitting 
transducers must be "looking " at the same area of the mean surface of 
the target for all relative angular positions to make the evaluation of 
aperture effect practically possible. The other system is the rotary table 
that is used to rotate the target through the intersecting beams of the 
transducers thereby providing the ensemble of surfaces necessary for 
the calculation of the average of the received power. The axes of 
rotation of the transducer positioning system and the rotary table are 
offset as shown so that the illuminated area of the target changes as the 
table rotates. 

The electronic system shown in block diagram form with the labora- 
tory parameters in figure III. 7 and photographically in figure HI. 8 is 
composed of three sub-systems: the transmitting sub-system, the 
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The Acoustic Tanks of the Remote Sensing Laboratory- 

Figure III. 2 



The Acoustic Tanks and Peripheral Equipment 
Figure III. 3 




Motor 



3 rpm 


Schematic Diagram of the Mechanical System 
Figure III. 4 



Transducer Positioning System 
Figure III. 5 
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receiving sub-system, and the data processing sub-system. The operation 
of each of these sub-systems in conjunction with the others is discussed 
in turn below. 

The transmitting sub- system consists of a pulse modulated sine 
wave power oscillator and the transducer which it drives. It is necessary 
to use a pulsed system instead of the simpler continuous wave system 
because the tank is not "dark. " The unwanted reflections from the walls 
of the tank and the surface of the water are gated out upon reception. 

These reflections, or reverberations, manifest themselves as a type of 
noise after several bounces; however, this was not a problem here because 
of the high signal level. The pulse width (PW) of the modulating pulse 
was chosen to provide the beam width limitation necessary to the measure- 
ment technique. The width of 600 microseconds gave a discernible 
steady state condition for the extreme case when the depression angle 
was $ j = 70° . The pulse repetition frequency (PRF) of the modulating 
pulse was chosen to produce a time series of statistically independent 
received pulses. This was a choice made relative to the choice of 
angular velocity of the target which was the low value of 3 rpm to minimize 
turbulence of the water. For this value of 3 rpm the PRF of 5 cps was 
sufficiently low to produce independent data points for the two targets 
used in the experiment. 

The source of acoustic energy was a piezoelectric crystal transducer 
of the piston type, shown in figure III. 9. The transducers used were supplied 
by the Branson Instruments Company and have a nominal diameter of 7/8 
inches. The antenna patterns of the transducers used (several transducer 
failures occurred) were measured and each was found to have the same 
pattern (shown in figure III. 10) within the limits of experimental error. Also, 
the mechanical and acoustical axes of each transducer were the same. 

The antenna pattern predicted by theory for a piston transducer is [Morse, 
19481 
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Rotary Table and Drive Mechanism 
Figure III . 6 



Block Diagram of the Electronic System 
Figure III . 7 
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Principle Components of the Electronic System 
Figure III. 8 



Piston Transducers 
Figure III. 9 







Block Diagram of the Data Processing Sub-System 
Figure III. 11 
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power averaged over the ensemble of surfaces generated by the rotating 
target as each pulse in the time series is energy which has been scattered 
from a part of the rough surface sensibly different from the remainder. 

The calculation of the average of the time series of RF signals in 
this way is achieved through a conversion of the series to a continuous 
waveform from which the time average of the power is obtained. The 
series of signals and the waveform to which they are converted are shown 
in figure III. 12 (note the effect of beam width limitation). 

The conversion process begins by detecting the RF signal and then . 
sampling each detected pulse during its steady state portion, thereby 
determining the voltage samples from which the average power is computed. 
The time of sampling is determined by the arrival of a gate pulse generated 
in the oscilloscope to intensify the screen during the chosen portion of the ' 
waveform. The portion chosen here is from that part of the signal which is in 
the steady state., The series of voltage samples is converted to a series of 
dc voltage levels with each voltage level being equal to that of the 
corresponding sample from which it is derived; This series of dc levels 
is then made the input to the analogue computer which computes the time 
average power of the waveform 

T 

— fv^dt: 

T J 

O 

which is equivalent to computing the arithmetic average of the power of 

each of the voltage -samples of dc levels. 

The practical application of this scheme of calculation was made 
difficult by the fact that the data was taken at a point physically removed 
from the analogue computer by approximately one mile. This separation 
made it necessary to telemeter. the series of dc voltage levels to the 
analogue computer and then telemeier back the resulting answer. The 
information was sent, over the telephone lines by means of a telemetry 
system designed for the general transmission of analogue data. 

The calculation of the average power was done on a relative basis 
as the system was normalized with respect to a calibration power level. 

This calibration power was the power received by the receiving transducer 
upon transmission of a given power level by the transmitting transducer 
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Conversion of RF Waveforms for Data Processing 
Figure III. 12 
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one meter away. The procedure used in taking data was to read all 
calculations of average power from fading (scattered) signals relative 
to the average power obtained with the non- fading, calibrating received 
signal as the input. By following this procedure, the gains of the 
transducers and the transmitted power level were eliminated from the 
calculation of a Q , and eq. (Eli. 4) was used directly with r n = 1 because 
it is <^P n {$ Q , $ ^ , <j) 2 ^ that was measured. 

The measurement of the average received power is clearly a 
complicated process and therefore one in which the likelihood of experi- 
mental error is high. The error occurring in the measurement process can 
be put into two independent categories: those generated by the electronic 
system and those which have their origin in the mechanical system. The 
experimental error due to the electronic system is believed to be negligible 
because of the many precautionary measures that were taken. Before 
and after each measurement run the receiving and data processing equip- 
ment were calibrated on the power level of ithe .power oscillator which 
was monitored throughout the experiment. If any change in calibration 
occurred, the data was discarded. The accuracy of the data processing 
equipment was checked several times by recording the fading signal on 
film and performing a manual calculation of the average power from which 
to make a comparison check. In each instance the agreement was good. 

The greatest source of experimental error is believed to be the mechanical 
system, and this was due to transducer pointing error. It is necessary 
that the transducers "look " at the same area on the mean surface (for 
the way aperture effect is evaluated) and this requirement is difficult to 
fulfill because of the small beamwidths. The transducer mounts were 
optically aligned, and the alignment was checked before and after each 
measurement run. The post-experiment alignment was frequently found 
to be poor with always some minor deviation which had an effect which 
is difficult to evaluate. 

To insure the quality of the experimental data, runs were repeated 
until agreement within a db between two runs was obtained. It was, 
however, not possible to do this for the data on which the results of figures 
III. 23 and III. 25 of Section III E are based because the experimental 
surface became badly cracked and unusable. These results are, therefore, 
possibly in error. Other results for this surface shown in figure III. 21 
are, however, based on repeated data and presumably good. 
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III. D Description of the Rough Surfaces 

The characteristics of the re-radiation of waves from surfaces is 
determined, in general, by the properties of two random processes. These 
processes are the variation of surface heights and the variation of material 
parameters. For most surfaces, natural and artificial, these processes 
are of such complexity they can be described in only a very approximate 
way if at all. The problem of description is complicated by the lack of 
statistical independence between these two processes and their non- 
stationarity in space. A further complication is introduced by the possibility 
of a variation in time. Consider, for example, the problem of describing 
farmland on a large scale. Here the lack of statistical independence is 
shown by the fact that vegetation growing in valleys' differs from that on 
the crowns of hills. The non-stationarity of both processes in space is 
evident from the variety of farm crops with different surface textures and 
material properties and the variation in time is caused by the changing 
weather conditions. Because little. is known about the nature of spatial 
non-stationarities, the effects of the weather, and the dependence of 
vegetation growth upon surface structure, it would not be possible to 
obtain a description of the terrain even if there were a way to treat the 
statistical features of vegetation. 

/ 

There are, however, natural surfaces that are more amenable to- 
mathematical description. The surface* of the ocean has received much 
attention [Pierson, I960] and progress has been made to the point that a 
prediction of radar return is now possible which takes into account, to a 
degree, the statistical features of the ocean's surface*. This case is 
much less difficult because of the simpler structure of the surface and 
also because the material variation is caused, primarily, by foam on the 
surface which is probably negligible for many problems. However, the 
spatial non-stationarity and the variation in time exist, and these are 
the phenomena which are of great interest to oceanographers. Therefore, 
even in the simpler cases, there is great difficulty in obtaining a description 
of the surface which has some degree of completeness. 


* Work now in progress at the Remote Sensing Laboratory of the 

Center for Research, Inc. 
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The surfaces used in the experiments were chosen not to have 
the general complexity of a surface generated by two random processes 
and the material variation was eliminated by using homogeneous materials. 
The goal was to construct stationary surfaces of roughnesses that could 
be described mathematically and also be smooth enough so that, at the 
frequency used, the locally flat approximation would be valid. To an 
acceptable degree, this was achieved. 

Two surfaces of different roughnesses and different reflectivities 
were constructed. One of these was made by striking mild steel sheet 
(of the kind used to repair automobile bodies) randomly with the ball (of 
approximately 3/ 4 inches) of a ball pean hammer and then gluing this to 
an aluminum sheet using a dense paste to eliminate air-pockets. The 
other surface was made by flowing grout over a sand surface which had 
been smoothed by coating it heavily with fiberglass. The sand surface 
was made by gluing sand particles to an aluminum sheet. The presence 
of the sand destroys, to a degree, the homogeneity of the surface 
material. The reflectivities of the surfaces were not measured, and 
because the targets were not infinitely deep, the reflectivities near 
normal incidence are difficult to predict. In fact, because of the type 
construction used in making the targets, the problem is not unlike the 
layered media p roble m with the surface layer being rough. 

The differences in the materials used to make the targets appear 
as differences in the critical angles for the two surfaces. The critical 
angle is the angle beyond which all energy is reflected and none refracted. 
At this angle the quantity 


i 

/A -S unjo- 
in eq. (II. 5) becomes zero and with a further increase in 9 becomes imaginary 
causing J_p { to become unity. Estimates of the critical angles for the 
two surfaces, were made by estimating the value of the phase velocities of 
the shear waves, 4 |x '/p ' , in the materials by using representative values 
of the phase velocities of the longitudinal waves, 4 (X , ) v +2p, 1 , and 

assuming the Poisson relation, jx 1 = (k 1 )' > to hold. P The representative 
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values for the phase velocities of the longitudinal waves were taken from 
the Handbook of Physics and Chemistry [1960-1961] and are 5000 
meters per second for steel and 4000 meters per second for the rock- 
like grout. Using a phase velocity of 1500 meters per second for 

water, the critical angles are 0 =30° (measured from the vertical) for 

c 

the steel and 0 =40° for the grout, 
c 

The properties of the surface roughnesses of the two targets were 
investigated by making estimates of the probability distribution function 
of surface heights, the autocorrelation function of surface heights as a 
function of lateral distance, and the stationarity of the processes. These 
estimates were computed from series of sampled height measurements 
that were taken along profiles of each target. The measurements were 
made using a depth gauge mounted on a machined surface and driven by 
a lead screw mechanism which allowed the sampling interval to be 
varied. This apparatus is shown mounted over a target surface in figure 
III. 13. 



Apparatus for Making a Profile of a Rough Surface 
Figure III. 13 
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The sampled height measurements were made along several 
profiles of each target. The profiles were separated a distance sufficient 
(several correlation distances) to insure independence of the data points 
of separate measurement runs. The probability distribution function and 
autocorrelation function were computed for each of these runs as well as 
the mean, m, and standard derivation of heights, a. Similarity of the 
computed results from different parts of the same surface then gives an 
indication of the stationarity of the process over the surface. A serious 
difficulty encountered in making an estimate of the stationarity and the 
other functions as well was caused by the presence of large scale irregular- 
ities created in the construction of the targets. These irregularities, which 
created what are called "regional slopes, " were of a lateral extent large 
relative to the variations of the surface and had the effect of causing the 
deformation of the mean plane into some other unknown surface. The 
"regional slopes" were not great enough to cause serious error in the 
measurement of the scattering cross section but did prevent, for some 
discarded profiles, an accurate determination of the small scale structure 
superposed on them. 

The estimate of the correlation function was calculated using the 
statistical estimator (the overbar indicates sample function) 





r(i£x)= 


i 



III. 5 


where: h. is the i** 1 height in the series of N points. 
1 is the lag . 


A x is the sampling interval. 


N 



i=l 
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— 2 

The denominator of eq. (III. 5) is the sample variance, a , or the square 


of the sample standard deviation. The maximum lag used in the calculation 
of r , as recommended by Blackman and Tukey £19583 * was N/10 where 
N was made sufficiently large for the particular sampling interval used 
to adequately define r . The criterion for the selection of the sampling 
interval, according to the sampling theorem [Bendat, 1958], is that it 
must not exceed half the period of the highest frequency present in the 
record of the profile. For the surfaces used in the experiment, the value 
of Ax = 0. 01 inches amply satisfied this criterion. The number of points 
taken in the measurement of both surfaces was 1000 which corresponds to 
a record length of ten inches. This provided a maximum lag of one inch 
which is sufficient to adequately define the falls of the sample autocorrelation 
coefficients. The one exception to this was a measurement made over the 
steel surface with a lag of 0. 1 inches and a record length of 17 inches. 

This measurement was made to investigate the behavior of the autocorrelation 
coefficient at longer lag distances. The sample autocorrelation coefficients 
computed (by machine calculation) are shown together with the sample 
variances for each of the surfaces in figures III. 14 and III. 15. 

The behaviors of the sample autocorrelation coefficients indicate 
by the smooth falls and following rises that the surfaces are gently 
undulating and, to a degree, periodic. The rises could not be thoroughly 
investigated in the case of the steel surface because of the regional 
slopes difficulty. The periodicity is evidence of the artificiality of the 
surfaces as the estimated period of the low frequency components was 
approximately the spacing between the surface variations. This estimate 
was very gross for the steel target but easily made for the grout surface. 

The presence of the periodicity indicates that the sample autocorrelation 
coefficients are probably not well described by any simple function C Pierson, 
1960] . However, since this low frequency behavior is dominated strongly 
by the falls of the autocorrelation coefficients from the origin it is this 
region of the curves near the origin that was investigated. It was found 
that near the origin the sample autocorrelation coefficients are closely 
approximated by the functJ 


although this is not shown in the figures. However, away from the origin, 
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all the curves but one fall off more rapidly than this and a better overall 
approximation is the Gaussian function 



The probability distribution functions of the sampled heights were 
determined for each of the profiles for which a sample autocorrelation 
function was calculated. However, no calculation was done for the 
profile of the steel target made with the larger sampling interval. These 
results are shown in figures III. 16 through III. 19 together with plots 
of the Gaussian distribution function 


J 

-CO 



2 

obtained using for m and cx values of the sample mean, m , and the 
2 

sample variance, a . Comparison of the experimental curves to those 
calculated shows the near-Gaussian nature of the random height process 
as sampled from the four profiles. 

The estimates made of the autocorrelation function and probability 
distribution function are very similar from profile to profile for both of the ex- 
perimental surfaces and the processes are considered to be, for practical pur- 
poses, stationary. Therefore, because of the Gaussian behavior of the sample 
distribution functions, the joint probability density function is, for 
both surfaces, approximately Middleton, 1960 


$ [h,(x v y,)j ha.O^yj] = 


i 

<3iT \j | - r 3 - 



where, because of stationarity , m i = m 2 = m ' o -^ = a 2 = a and r depends 
only upon the relative distance between (x^ , y^) and (x 2 » y 2 )* The values 
of a , m, and r used in this expression are averages of the values obtained 
from the measurements of each profile. 
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Since the random height processes for the two surfaces are, for 

9 h 

the purposes at hand, Gaussian and stationary, the surface slopes, -g-y 
where 1 is the distance along a profile of the surface, are, according 
to Middleton [ 19603 , normally distributed with mean zero and variance 


— <~j T 


x =o 


For a Gaussian autocorrelation function, the variance is 


tf 9 - 

l*- 

Although, the sample autocorrelation coefficients are clearly not fit 
well by the Gaussian function, especially near the origin, the variance 
of slopes can be estimated by using the values of sample variance and 
the calculated L’s, as given on figures III. 14 and III. 15. For the steel 
surface, the values of the variance of slopes for the two profiles made 
with the shorter sampling interval are estimated to be 

4.89 xlO -3 3.72 xlO -3 

with the average of these being 4. 30 x 10 . For the grout surface 

this estimate is more difficult to make because one of the curves is 
closely fit by the three halves power function. The problem is circumvented 
by sending a Gaussian function through the point r = l/e , and using the 
corresponding value of L (0. 26 inches) in the calculation of the variance 
of slopes. Doing this and proceeding as above for the other experimental 
curve gives the results 

12. 16 x 10 -3 6. 84 x 10 -3 

with the average being 9.5x10 . Therefore, on the average, the grout 

surface is approximately twice as rough as the steel surface using the 
variance of slopes as a criterion. 
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HI . E The Experimental Results and Comparison with Theory 

The experiment to determine the average differential scattering 
cross section as a function of the angles (<t> Q > <t > } / <t> 2 } was performed in 
a way designed to show the characteristics of two distinct features of 
rough surface scattering. One of these is the nature of the scattering 
when the transmitter and receiver depression angles are held equal 
throughout their variation (</> .. = <$> ^). This includes scatter in the specular 
direction ($ Q = 0° ) and also backscatter 180° ). The measurement 

of scatter in the specular direction as a function of depression angle is 
important because the differences in surface properties, both material 
and roughness, are so strongly evident there. Backscatter is, of course, 
important for practical reasons and also because of the large amount of 
accumulated data with which comparison can be made. The other feature 
of omnidirectional scattering which influenced design of the experiments 
is the nature of the scattering for one transducer held fixed in angular 
position (<£ 1 = constant) and the other variable. The behavior in the plane 
of incidence (<t> Q = 0° ) under this condition is of interest because a 
comparison can be made between the results of La Casce and Tamarkin 
Cl 95 61 for scattering from regularly rough surfaces and those obtained here 
for randomly rough surfaces. 

The results of the measurements made with the receiver and 
transmitter depression angles equal are shown in figure III. 20 for the 
steel surface and figure III. 21 for the grout surface. These are curves 
with a Q plotted against with $ } ? varying parametrically. Comparison 

of the results for the two surfaces shows, at least on an intuitive basis, 
that the grout surface is rougher than the steel. This is seen from the 
smaller values of cr c measured in the specular direction and also from the 
slower decrease of a Q with depression angle for angles away from the 
specular direction, although this is not as evident. The values in the 
specular direction for both surfaces are seen to be nearly constant with 
no discernible pattern to the small (on the db scale) changes. The results 
for the steel surface have an anomaly that is not present in the grout 
surface results. This is the crossing of the curves for angles <j) 1 = <f> 2 = 

50° , 60° by the curves for iower depression angle <t> t - 6 0 = 30°, 40°. 
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This is not believed to be the result of an experimental error because the 
experiments were repeated several times to observe this phenomenom. 

The experimental results for $ ^ held constant at 40° and <f> ^ 
variable are shown in figures III. 22 and III. 23 for the steel and grout 
surfaces, respectively. For these curves, cr Q is again plotted against 
<t> 0 with <j) 2 as the parameter. To show some features of these results 
with greater clarity cuts of 4> D ~ constant were made in figures III. 22 
and III. 23 to obtain plots in which a Q is plotted against </> ^ with as 
the parameter. These curves are shown in figure III. 24 for the steel 
surface and figure III. 25 for the grout surface. The results given in 
figures III. 22 and III. 24 for the steel surface are seen in greater per- 
spective in the three dimensional representation of these results shown 
photographically in figure III. 2 6. The three dimensional figure is 
formed by the radius vectors emanating from the intersection of the 
principal rays of the antennas, which is the origin of the coordinate 
system of figure III. 1. The length of a radius vector is proportional to 
the a Q in db for the angles of observation ( <f> Q , <j) in which the vector 
is oriented. The principal ray of the transmitting antenna which defines 
the angle of incidence $ ^ = 40° relative to the normal of the mean plane 
of the rough surface is depicted by the shaft piercing the figure. The 
view of figure HI. 26a is that seen from the sidescatter position, <f> ^ =90°. 

The views of III. 26b and III. 26c are those seen from 6 =45° and 

^ o 

= 135° , respectively. The results for the grout surface are not shown 
three dimensionally because of their similarity in form to those for the 
steel surface. The principal difference between the results for the two 
surfaces is that the peak of a Q in the specular direction is more pronounced 
for the steel surface which indicates that the grout surface is the rougher 
of the two. The relative roughness is shown in another way by the 
vertical extent of the curves of figures III. 22 through III. 25; for the steel 
surface the values of a Q vary over a 35 db range while the variation for 
the grout is only about 20 db. A peculiar feature which distinguishes the 
curves of figure III. 24 from those of figure III. 25 is the behavior far 
away from specular U Q > 60° ). For the steel surface there is an increase 
and then slight decrease as the depression angle increases while the 
grout surface has a general decrease with depression angle with a hint 
of an increase near <f> 2~ 30° . It is recalled that these results for the 
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grout surface are suspect (see Section III C) and perhaps no great amount 
of credence should be placed in this peculiar difference in the curves. 
However, considering the anomalous behavior of the steel target results 
mentioned above it is not unlikely that this is a similar phenomenom. 

The experimental results shown in figures III. 20 through III. 25 
were computed using values of re-radiated power that were assumed to 
be predominantly scatter. It is recalled that this is, in fact, a condition 
for computing a differential scattering cross section. The degree to which 
the power was being reflected was not determined and a separation into 
scatter and specular components is not possible. However, from Section 
II C it is seen that the factor of the specular power density 


ex 




III. 6 


can be used to indicate to what degree reflection is occurring. Using 
average values of the variances given for the two surfaces in figures 
III. 14 and III. 15, the factor III. 6 becomes 



III. 7a 


for the steel surface and 


3 . I (cos cos 

for the grout surface, where the conversion <j> 2 = e oz * $ i = 180° - 0^ z 
has been made. For the depression angles $ j = ^ 2 = ^ it is seen that the 
quantities III. 7 are quite small. However, for larger depression angles 
they increase in size and at<^ = $2 = 70° (the limit of the range of 
angles over which data was taken) become 

-4. 23 
e 

for the steel surface and 

-1.44 

e 

for the grout surface. The factor for the steel surface is still quite small; 


III. 7b 
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Figure III. 20 
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EXPERIMENTAL AVERAGE DIFFERENTIAL SCATTERING CROSS SECTION FOR THE GROUT SURFACE 
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experimental average differential scattering cross section for the steel surface 

VERSUS AZIMUTHAL ANGLE *0 FOR ♦i=40°WITH *2 AS THE PARAMETER 
Figure III. 22 
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experimental average differential scattering cross section for the grout surface 
VERSUS AZIMUTHAL ANGLE ♦<> FOR 4| = 40° WITH*; AS THE PARAMETER 
Figure III. 23 
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VERSUS POLAR ANGLE <t> 2 FOR'Ci = 40°WITH$iMS THE PARAMETER 

Figure ill. 24 
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EXPERIMENTAL AVERAGE DIFFERENTIAL SCATTERING CROSS SECTION FOR THE GROUT SURFACE 
VERSUS POLAR ANGLE ! *2 FORI <h = 40 ° WITH 1*0 .AS THE PARAMETER 


Figure III. 26a 



Figure III. 26b 



Figure III. 26c 


Three Dimensional Representation of Omnidirectional Scattering 
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III. 8 

From this it is seen that the range of bracketing must extend from the 
origin where r = 1 to a point where r is small enough to make the 
quantity negligible. The values of a used in III. 8 in the determination 
of the range of r and those calculated for III. 6. Substituting these 
gives 


<2Xp 


III. 9a 


for the steel surface and 


€Xp- -3.1 0~ r )j- 


III. 9b 


for the grout. It is seen that the range of bracketing depends upon the 
angles of incidence and observation. To insure bracketing over the range 
of angles for which there is experimental data, the angles for the steel 
surface are chosen to be <^ = $2=70°, and III. 9a becomes 


exp 


£-4. 23(l-r)| 


III. 1 Oa 


A similar choice cannot be made for the grout surface because the 
factor of (1 — r) becomes too small (from earlier, it is 1. 44) for even r = 0 
to make the exponential negligible, as required. It is this inability of 

t 

the exponential to become negligible that leads to the appreciable specular 
component at these angles, as discussed earlier in this section and in 
Section II D. A marginal choice for the grout surface is $ ^ = $ 2 = 60° ; 
for this set of angles III. 9b then becomes 


exp 


£-3.1(l-r>y 


III. 10b 
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For the exponents of III. 10, the range of bracketing is chosen to be from 
r = 1 to r = I/e; this choice is motivated in part by the fact that the 
Gaussian functions of figures III. 14 and HI. 15 pass through these points. 
Therefore, the correlation distances for both the exponential and Gaussian 
function are those lag distances corresponding to the points r = l/e on 
the curves of sample autocorrelation function; for the two profiles made for 
each surface, these values are 
L = 0. 26, 0. 23 inches 
for the grout surface and 

L = 0, 55, 0. 70 inches 
for the steel. 

The other difficulty encountered in applying the theory is that 
the reflection coefficient is unknown near normal incidence, as mentioned 
above. Away from the vertical, past the critical angle, the reflection 
coefficient approaches unity and as it does so the surface becomes 
perfectly rigid. Because of this, the results of Chapter II are used under 
the condition that the surface is perfectly reflecting, and comparison is 
made between theory and experiment for angles past the critical. 

The expressions used in the comparison are eq. 's (II. 28) and (11.29) 
after specialization to the condition of a perfectly rigid surface. This 
specialization is made through the a's which become for X 1 = 1 

a. - 2k cos 9 
1 oz 

a n = -2k cos0 

u OX 

a 0 ^~2k cos0 
3 oy 

Substituting these values and the expressions for the K's which are 
K x =k(oose ox -cose lk ) 

K y «k(cose oy -cose ly ) 

K z = kfcose oz - C°s6 lz ) 

into eq. 's (II. 28) and (II. 29) and manipulating gives the results 
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hi. 11 


for the Gaussian autocorrelation coefficient and 
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cr 0 - Jcos^o^ + 9 (iios^ -cos^x) cos e ox co^©ox' 

V. (cos - 0 OZ. — COS &12_J 

-+ a. cos © oz coj ©by 
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III. 12 

for the exponential autocorrelation coefficient. The conversion from the 

direction cosines of the directional angles (© ,0 ,0 , ©. , 0. , 0. 

* ' ox' oy' oz' lx' ly' 1 

to the spherical system ( (j> Q , <f>^, $ is given by 


cos©., = 0 

X X 

cos© = c osd> o 
oz v 2 

cos© jy = sin</> 1 _ 

cos© Q y = cos(J» o sin<^ 2 

COS©-. = ~COS(f) . 
lz i 

cos© ox = sin$ o sin^2 


The values for the variance of slopes used in the expressions for a Q for 
the Gaussian autocorrelation function are those average values determined 
in Section III D for the two surfaces. The ct q for the exponential auto- 
correlation function depends upon the quantity 
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r 


L/ ?\ 


i 3-~ 


L<°*rJ 


The values of this quantity used in eq. III. 12 are averages obtained using 
the variances for each surface given in figures III. 14 and III. 15 and the 
correlation distances given above. These average values are 

2.54 xlO 3 

for the grout surface and 

1.872 x 10 3 
for the steel surface. 

The theoretical results for the exponential correlation function are 
shown, together with the experimental curves, in figure III. 27 for the 
steel surface and figure III. 28 for the grout surface for the' angles j " 

2 = *0° / 30° , 50° , 70° . For the angles <j> ^ ^ ~ 20° * 40° , 60 6 

the results are shown in figures III. 29 and III. 30 for the steel and grout 
surfaces, respectively. A comparison between the experimental and 
predicted values for the steel surface shows that there is closer agreement 
for the angles $ ^ = $ 2 = 50° , 60° , 70° with the best agreement occurring 
for the two larger depression angles. For all smaller depression angles the 
agreement is not close especially in the neighborhood of the specular 
direction. As backscatter is approached, the agreement improves with 
the only serious difference occurring for the anomalous curves for the 
angles $ i ~ <!> 2 = 30° > 4 0 0 . A comparison for the grout surface shows 
fair agreement between the predicted and experimental results but near 
specular there are wide differences especially for the larger depression 
angles where agreement should occur. It is clear that for both surfaces 
the theoretical results obtained using the exponential autocorrelation 
function do not generally give accurate predictions of the experimental 
results for angles in and near the specular direction. This is not surprising 
as the sample autocorrelation functions are not fit closely at all by the 
exponential function. What is surprising is the closeness of the agreement 
between theory and experiment away from the specular direction near 
backscatter. 
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The predictions of the experimental results made using the 
Gaussian autocorrelation function are shown in figure III. 31 for the steel 
surface and figure III. 32 for the grout for the angles = , 30° , 

50° | 70 ° .• The results for the angles <f> ^ ^ ~ 20° , 40°, 60° are 

shown in figures III. 33 and III. 34 for the steel and grout surfaces / respectively. 
The predictions for each of these surfaces are shifted downward; for 
the steel surface each curve is smaller by the factor 0.5 than that calculated 
and for - the grout surface the factor is 0. 22. It is seen that by doing this 
close agreement is obtained for both surfaces for angles that do not differ 
greatly (about 25° ) from the specular. However, far away from the 
specular direction, the theory disagrees greatly with the experimental 
results and is of little value. For angles less than the critical angle it 
is expected that the predicted values for X 1 = 1 would have to be shifted 
downward. However, for angles beyond the critical this should not be 
necessary and there should be agreement on an absolute basis. Since 
the agreement is only relatively close, the use of the Gaussian auto- 
correlation function is seen to lead to error in the specular direction as 
well as away from it, although to a much lesser degree. Also, the error 
is greater by about the factor two for the results for the grout surface 
than for those of the steel surface, which is expected as the steel surface 
is, from Section III D, described more closely by the Gaussian auto- 
correlation function than is the grout surface. 

The theoretical results for $ l constant at 40° are shown in figures 
III. 35 and III. 36 for the steel and grout surfaces, respectively, for the 
exponential autocorrelation function. The general agreement for both 
surfaces is clearly not good. The results for the steel surface give rather 
close agreement for the angles 2~ ' ^0° ' e l sew here there 

are wide differences. It is recalled that the agreement in figures III. 27 
and III. 29 is also closer for these angles than for others. The theoretical 
results for the grout surface would be closer to the experimental values 
if they were to be shifted downward by 5 db. A shift in either direction 
would bring no improvement to the results for the steel surface. It is 
important to note that the theory predicts for the grout surface the general 
decrease of cr with increasing depression angle, for angles far 
away from the specular direction, as discussed earlier in connection with 



87 


figures III. 24 and III. 25. The increase of a with depression angle 
observed for the steel surface is not accounted for by the theory. 

The predictions made using the Gaussian autocorrelation 
coefficient are shown in figures III. 37 and III. 38 for the steel and grout 
surfaces, respectively, for a constant of 40° . The values for both 
surfaces are shifted downward by the same factors used for the results 
shown in figures III. 31 through III. 35. By doing this, fairly close 
agreement is found for angles that are near the specular but elsewhere 
there are great differences; this is the same result as found earlier for 
the Gaussian function. The values for the steel surface compare more 
favorably with experiment than those for the grout. The predictions for 
the grout surface would be in closer agreement if all values were shifted 
downward an additional 3 db; that this is so indicates the possibility of 
experimental error, which is suspected for these results. 

In the comparison of the experimental results shown in figures 
III. 20 through III. 25 to the predicted values for both the exponential 
and Gaussian functions, several important features of the agreement 
between theory and experiment become apparent. The results based on 
the exponential autocorrelation function are seen to agree with experiment 
over a wider range than do those for the Gaussian function. It is 
surprising that this is so considering how well the curves of sample 
autocorrelation coefficient are fit by the Gaussian function. It'is true, 
however, that the Gaus sian-based predictions give a better description 
of the experimental results in and near the specular direction, especially 
for the steel surface which has a sample autocorrelation coefficient more 
closely fit by the Gaussian function than the grout. It is this goodness 
of fit in and near the specular direction which is the principal success 
in the use of this function; its most disappointing failure is its inability 
to predict the experimental results at other angles of observation. It 
is in these directions, those away from the specular, that the results 
for the exponential autocorrelation function are most successful. These 
agreements (and disagreements) of theory with experiment are seen more 
clearly in figures III. 39 and III. 40 in which the theoretical curves for 
both autocorrelation functions are plotted against the experimental results 
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for both surfaces for <£ i ~ $ 2 = ; resu ^ ts f° r Gaussian function 

are shifted downward as before. Comparison is made for the steel surface 
in figure III. 39 and the results for the grout surface are compared in 
figure III. 40. The results for the two autocorrelation functions are seen 
to bracket the experimental values for the full range of azimuthal angles 
with the exception of the results for the grout surface in and near the 
specular direction for which the experimental values fall below the 
theoretical ones. This bracketing, which occurs for other polar angles, 
gives strong evidence from the viewpoint of comparing theory with experi- 
ment that the true autocorrelation coefficient lies between the exponential 
and Gaussian functions. Since this is the case, it is concluded that 
the theory as based upon the Kirchhoff approximation is a valid one. 
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CHAPTER IV 

APPLICATION OF THE ACOUSTIC RESULTS 
TO ELECTROMAGNETIC SCATTERING 


Many of the radiation phenomena associated with different physical 
phenomena possess a high degree of similarity. A phenomenon such as 
diffraction is common to fields as diverse as acoustics and quantum 
mechanics. Advantage can be and has been taken of such similarities 
if they are sufficiently alike to permit analogue study. It is not even 
necessary for this to be true to benefit , for if theories for two (or more) 
physical processes have a common foundation then a verification of the 
theory for one process implies general applicability. The principal 
reason for performing model experiments whether for the verification of 
theories or for making analogue studies is economy of effort. 

Acoustic model experiments have long been used to study a variety 
of electromagnetic problems. Jordan £1941 3 ; in an early application, 
made model studies of the properties of antenna arrays and Maestri Cl 961) 
has modeled antenna radiation patterns. Meyer Q958J made use of the 
technique to study microwave surface reflection losses and EdisonCl963] 
has used it to investigate scattering in turbulent media. Recently, 

Ayiku and Moore Cl 9653 have performed model experiments on knife edge 
diffraction in which both single and double edges were used. Probably 
the greatest amount of activity has been in the area of simulation of 
radar backscatter, the reason for this being the difficulty and expense 
incurred in performing full scale measurements C Moore, 19623. 'Much 
of this work has been and is being done at the University of New Mexico 
where a number of simulation studies have been carried out. These 
studies include simulations of radar back scattering from several types 
of terrain including citiesC Warner et al. , 19 623 C Edison, 19613 , a 
simulation of the lunar radar echo£Hayre, 19623 and the modeling of 
the radar cross section measurement of geometric shapes CKoepsel and 
Ahmed, 19623 . A complete description of the modeling techniques used 
for radar backscatter as well as other electromagnetic problems is 
contained in the dissertation of Edisonf 19613 . 
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A major difficulty in making acoustical analogue studies of some 
electromagnetic problems has been that the limits of applicability were 
not clearly defined. The reason for this is that the extent to which the 
scalar acoustic wave simulates the vector electromagnetic wave had not 
been fully described. In some- problems the vector nature of the electro- 
magnetic wave is not important and can be neglected, a case in point 
being scattering in a turbulent medium [Edison, 1961, 1963] . However, 
in other problems this is not possible. Consider the reflection of a 
plane electromagnetic wave of arbitrary polarization from a plane surface. 
The reflection is a vector phenomenon as the two polarization components 
are reflected unequally [Stratton, 19411 , A simulation can be made of 
the reflection of the separate components but it is not possible to model 
directly the generally oriented vector reflection. 

The problem of the applicability of the analogue scalar measurements 
does not arise when the model experiments are being performed for the 
purpose of verifying electromagnetic theories that are based on assumptions 
identical to those of the analogous acoustic theory. This application of 
acoustic model experimentation is especially important for the study of 
the re-radiation of electromagnetic waves from rough surfaces as all 
electromagnetic theories have bases that allow the development of the 
simpler acoustic form. 

The applicability of the results of model experiments of acoustic 
wave re-radiation from homogeneous, rough surfaces to the analogous 
electromagnetic problem is determined through the use of theory based on 
the Kirchhoff approximation. Of course, this determination is accurate 
only for situations for which this approximation applies. The limits of 
applicability are determined by comparing the expressions for the acoustic 
re-radiation field with the expressions for the analogous electromagnetic 
re-radiation field and determining under what conditions the vector re- 
radiation differs enough from the scalar to destroy the similarity. The 
acoustic expressions used in the comparison are those developed in 
Chapter II; the electromagnetic expressions are developed below. 

The derivation of the expressions describing the re-radiation of a 
plane electromagnetic wave from a rough, homogeneous surface proceeds 
similarly to that for the acoustic case, the principal difference being the 
complexity caused by the vector re-radiation at the rough surface. A 
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form of the vector Helmholtz integral is used with values of the electric 
and magnetic fields tangent to the surface being determined through the 
use of the Kirchhoff approximation. This problem has been treated in this 
way before and the presentation here is similar to that of Kovalev and 
Pozdnyak [ 19613 , the difference here being the inclusion of a two 
dimensional rough surface. 

The re-radiated electric field, E (P), at the far field observation 
point P (see figure IV. 1) is given by the form of the Helmholtz integral 


E(P)~ ''“'Mo* 

■4-~rr R 0 




kR 0 



TtxH -n (Tax Tt. rV) 



IV. 1 

where r Q is the unit vector in the R q direction, |i ^ and e are the per- 
meability and permittivity of the medium, respectively, S is an aperture 
of the size discussed in Section II. B and n is the unit normal to the 
surface. This form' is derived in Appendix V; a derivation apparently 
appears in the Russian literature [ Frandm, 1957] but this was not available. 
For the purpose at hand, the factor of r^ in the second term of the inte- 
grand is neglected as components of the fields aligned in this direction 
of propagation carry no power in this direction. 

The Kirchhoff approximation is used to determine the values of 
n x E and n x H at each point of the rough surface z =^(x,y). This 
is done by resolving the incident plane wave into local polarization 
components parallel and normal to a local plane of incidence and treating 
the reflection of the two components separately. The incident electric 
and magnetic fields are given by 
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hi- = r| A^ X 06 <g 

f 7 ' 4 , 

where a is the polarization vector and rt. is the unit vector in the 
direction of propagation," h = n.k . The resolution of and into 
local polarization components is facilitated by the establishment at 
each point of the surface a local coordinate system based on the plane 
of incidence formed by in and n . The unit vectors of this system are 
(see figure IV. 2) 

"n 

h = '0 o *n‘ 

5in-e IV. 3 


where © is the local angle of incidence and a function of (x,y) through 
n , as described in Section II. B. Resolving the incident field into 
local polarization components using the local coordinate system gives 
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for the parallel polarized components. 



• The Local Coordinate System 
Figure IV. 2 


With the resolution into polarization components, the tangential fields, 
n x H and nXE , are determined following Stratton [ 1941} . Beginning 
with the normally polarized component, the expression for the total 
tangential electric field for this component at the surface is 

nxE x =lAx e JlO + -Cl) iv. 5; 

where the reflection coefficient is 



and the subscripts 2 denote the parameters for the homogeneous medium 
upon which (E^, FL) is incident. The total tangential magnetic field for 
the normally polarized component is given by the following calculation 

"n x Hj_-= x + 1 a x Hj[ Ou* ^ (n r x. ^ 
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V^'4, 


IV. 5b 


where (Ej, , Hj_ ) is the reflected field, n is the reflected wave vector 
and 

rt. • n = -n • n 
1 r 


The parallel polarized components are treated similarly. The expression 
for the total tangential magnetic field is for this polarization component 



IV. 5c 


where 
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The total tangential electric field for this component is given by 

rtx £^+ nx E// =^x ( H//' +fp< ( H/j x V0p)^ | /£ l 
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The total tangential electric field is given by the sum of the tangential 
fields of the two polarization components, or 

nxE=nxE x + nx-E^ ■ IV. 6a 


and similarly the total magnetic field is 
rf x H = rt x H x + nxl^, 


IV. 6b 


where the right hand sides are given by eq. (IV. 5). 

The solution for the electric field at the observation point P is 
obtained by substituting eq. (IV. 5) into eq. (IV. 6) and substituting this 
result into eq. (IV. 1). The expression which results from these sub- 
stitutions is quite lengthy for the general case and is not given. Instead, 
expressions are developed- for special cases under several conditions of 
approximation for comparison with the analogous acoustical expression. 
These cases are those of normal and parallel polarization relative to the 
plane of incidence defined by the mean, normal, Ic , and FT . The acoustical 
quantity with which comparison is made is the pressure and the expression 
used is obtained from eq. (II. 9)- by neglecting second order partial 
derivatives of the surface and writing the reflection coefficient in closed, 
form. This is 



-f A; ('feo'lA ) ( H - -P ) <2- d S 


IV. 7 

It is recalled that eq; (II, 9) was obtained under the condition of a plane 
wave incident upon a rough aperture of homogeneous material. 

The acoustic quantity used in the comparison is chosen to be the 
pressure for reasons of convenience. The particle velocity can, of course, 
be used; however, by judiciously choosing the analogous electromagnetic 
quantity, this calculation can be spared. Edison [l96l] has shown that 
for the case of parallel polarization the tangential component of the incident 
electric field and the normal component of the incident particle velocity 



. Ill 


both vary as cos0 and, therefore, are necessarily analogues. The 
analogy is complete for this case because the incident magnetic field 
and the incident pressure do not vary with 0 (except in the phase factor 
which is common to all quantities). For normal polarization the situation 
is reversed and the pressure is analogous to the electric field. Therefore, 
for the case of parallel polarization the expression for the re-radiated 
pressure field is compared with that for the re-radiated magnetic field 
(to be developed subsequently) and for normal polarization with the 
electric field. 

The incident field for polarization normal to the y-z plane is 


~'aS~ 




Mu 


where n^ = j cos0^ + k cos0^ z . For rT. in the y-z plane the unit vectors 
are, at a point (x,y) on the surface 
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as before. These expressions are clearly complicated; to obtain a simple 
first order approximation, Z x and Z^. are set equal to zero (smooth surface). 
Then 

f=T p=7 n=lc 
as expected and, therefore, from eq. 's (IV. 4) 

— i " v'C' "fi 

E;£ = ;£. 

til = 

A 'r r 
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From eq. 's (IV. 5) and (IV. 6) 

fixt = 0 + -QJ 

= fix* E 

IV. 9a 


and 
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. 4=* -i 

r^<r p 


^'4 


IV. 9b 




Substituting eq. *s (IV. 9). into eq. (IV. 1) and neglecting the factor of 
gives for the electric field at P 



-JL £e 




4tT £ c 


’//{•**. 0 ~£±) 
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+ (^•£o)£ (l +£^) ^ (\ 




IV. 10 
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where, as before, K = lc 0 - ^ and kT^ = iTk , kT Q = r^k . Comparison of 
this expression with the analogous acoustic expression given by eq. (IV. 7) 
shows that when the observation point P lies in the plane of incidence 
the integrals are identical since I* • — 0 . For this case, the electro- 

magnetic re-radiation is essentially scalar; there is no depolarization of 
the incident wave and it is necessary to give only the single reflection 
coefficient £j_ • Out of the plane of incidence, there is a change of 
polarization which is a vector effect that cannot be simulated by scalar 
means. However, it remains to investigate the possibility of simulating 
the amplitude of the vector field. To do this it is necessary to examine 
only those components of E(P) normal to r Q , the direction of propagation. 
These components vary as the sine of the angle between ~r Q and E(P). For 
the first two terms of the integrand of eq. (IV. 10), there is a difference 
in the variation of sin© from the similar terms in eq. (IV. 7). For the 
last term in the integrand of eq. (IV. 10) for which there is no comparable 
acoustical term, the power producing components vary as the sines. of 
their respective angles withT Q times T* 1c , which causes the disappearance 
of this term in the plane of incidence. Therefore, the exactness of the 

simulation of the re-radiated field becomes less as the vector r* moves 

o 

away from the plane of incidence whether polarization is considered or 
not. The degree to which this occurs is determined by the variation of 

sine ox andr -‘ 6 o =oose ox • 

The expression for the electric field given by eq. (IV. 10) was 

developed under the approximation that the fields were re-radiated as if . 

the surface were smooth. The effect of the surface roughness on the 

vector re-radiation is brought out through consideration of the simpler case 

of re-radiation in the plane of incidence ( cos@ n = cos© = 0). This is 

lx ox 

done under the less severe restriction that Z x , Z y <<1. The unit vectors 
are approximately from eq. (IV.‘8) 


n = k 
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The resolution of the normally polarized field into local polarization 
components gives, from eq. (IV. 4) 
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(~ZxZy COS Qty j- Zx COS 012_) jZ. g ‘ 

\f^ |^-j- 2)<-f- 4Ly~— (“ co<s 0ty + cos 0-i-z.} 

X "'^ c °' 3 &y ^yQ-VS&y-rZyCoS ©l2-)-i-Zx~ COS 0 ( -2_ 

— CO^ 0IZ. ^Zjf^CO.S ^ty + COS 0ry + ZyCO^Oi^jv 

[^ S ©iz.^-^’yZxco.s ©y 4- Z x cos <3h^)j 

+1c ■— COS 0y Q-ZyZ.'K COS Oy + Z X COJ -0- t2 _ ^ 


Substituting into eq. *s (IV. 6) and (IV. 5) gives 
Taxt ='V\x E^-v-nx £// 


. Ad (\ + 0 *) 

jj A-Zx>Zy~- (- 2y cos eiy V Cos ^z-) 3 ”. 

/CZx cos 0 - l2 _ -Kp Zy cos ©y cos 0\-2- 
-v- (cos aiz-)C^xZy COS 0<y 4-Zx cos 0-iz.yi ~ -0/) Erg ** \ P 
\ + Z^+Z^h- (- Zy cos Oy •+■ Cos ©1-2.) ^ 

6^-iy ^y + ^y 005 ^ lz ^) + Z^COJ 0-1 yj 

— CoS 012 . jzjcos 0 y -V Cos &\y 4 -ZyC oS • 


A-(Co.S 


+ 


-Zi 

J- 


COS 0"i2_ 


ZyZyCoS ©y 


CoS 0-<yHr 2 X CoS 


■i] 
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nx H —1 a x 


-Xt,. 




— y-^v co. s &\y + cg?.s €hz. ^ &<g 

V^, 2/ +■ Zy~— (— Zy'CO.S <&iy + C*j.S -©vz. 

A ^-Cose^(-^'Zy<cos e ^+2xc^jd fz _)^H-'^-|co^€>,y 

f/ ^ e / + Z / °°^ ^0+ OJ^l zj-CO-S e,z jz^CO^ ©m, 

’+ ws 0y i- 2^, co«s ' 


-h C-COS Olz. ^(cos Q-y + 2y COJ- e 12 .) ( l -_Q j £ e ~^ ‘^ 

^F%I J^+2x4-Z^~~- (- 2y cos O-y -+ cos 

^ (cos Byj-ZycoS 0tz.J^( z X^e,^ -f-~£(Z X COsOy) 

Substituting into eg. (IV. 1) and neglecting the factor of as before gives 
the result . r 

t( P)= ik£ 2 ^ ( r 

^TTRo J s 


4 cot ©iz. (coiOry + Zy COS ©-le^Ci -J^) 


-oos^iz^ZxZ^cuxse^ + z^oo&etz.)^^^] _ coso^cox©^, cot By 
(-2 x Zy ce?-s 0iy + 2* a^s a t -2L )®X’l ~J } /y J - coiteiz, cos 002. (~Z y co* ©iy 
~*-Zx cos a, z)^(| - ^ - W ©02. (cos /Lycos. Qx'z.'f'i} 

^U,s&^fXc a s&y + Z /W & lz )Z x (, ^ ^ 

-+ 2 X cos ©I 2 ^cote^ + 2^ co-s ©jz^Z^os ©hy £l +jJ,J 

+ C'ZxZ^ cos 0 y -fz* cor* &1JL f .jnj 

- (- ^ ey + C« s e, e^+ujie^, t2/WjSli j 

e,i + ' c « ©>*• ©0^ (-2x 2y CO* Sy + 2x co^ e, t ) 
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+ cos COS ©oz. cos (-Z * zy c os. Q-y + Zx COiSii) Zk~ 0-£//) 

~ cos^cosOo* (-Z xZ/ ^y + z „ c -s©,^(z^ d / " 

&.*) ("l -■£//) ■+■ (COS ©i^. + Zy CUJ -©12.) 

(COS Oos. Cos O l3 .) 2x (I fi |T_ )"!. 

4-t|(-COS«, i )(C0i0r X + ZyCOS ©, 3 .) 2y.COS (| -_Pj 
4-(uDvS©^+-Zy c^s-©i2.^-cosOoy cos -©,^) Zx 0 + -£j_) 

~ Co-S -O-IZ. CO^ Ooy OOS ©ly (—Zy.'Z.y COS ©I 2^ + Z* COS Oi2L^ 

{^6o>S ©y +- 2y Co ■S“Oi2_ ) Zy ^|— X]y ^ 

-cos©-,^ c^e oyco^ ey (-z^ZyaosBy -*- z X u,sQ t z) 

C^x COS e t 2.)0 ~J}y) + COS ©12. CO.X©oy cos ©12. (- Z*Zy COS ©iy 

H-Zx<^S©( 2.^ 4- (ZxcOS-©jy-f CO-5 &y +■ Zy CO j ©I2_^| ‘-J^j^.l 


XK-r 


l + ZxS- 2y~~ — (— ~Zy COS ©ty +■ cos: ©I z 3^ 


IV. 12 


The comparison of this result, eq. (IV. 12), with eq. (IV. 10) indicates 
the severity of the approximation made in developing this latter equation, 
thereby illustrating the effects of the vector re-radiation. One of these 
is the depolarization of the incident radiation by the rough surface; this 
is seen by the presence of the non-zero factors ofT and Tt, which are the 
cross polarized components relative to the polarization vector T of the 
incident wave. However, the depolarization components all have Z or 

X 
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Zy or both as factors which indicates their smallness relative to the 

polarized component which has terms that do not depend upon the slope 

of the surface. These slope independent terms also depend only upon 

and not upon •R/ , as expected. For backscatter at normal incidence 

(cos©, = -1, cos© = 1, cos© , = cos©„ - 0) eq. (IV. 12) reduces 

j- ^ oz 1 y oy 

(considerably) to 




£ 0 








IV. 13 


At nonnal incidence, the reflection coefficients have the same magnitude 
but differ in sign, and their slopes with respect to angle of incidence are 
zero [ Stratton , 19413 ; therefore, the cross polarized component disappears. 
The same result is obtained for the case of perfect conductivity of the 
rough surface for which £^= 1, _Q_= • 

The field re-rediated from a perfect conductor is obtained more 
generally by using the boundary conditions 

n x E = 0 

n x H = 2n x i. IV. 14 

which are obtained from eq. 's (IV. 5) using .Q, = 1 , _Q_= -1. For the 

case of normal polarization, substitution of eq. (IV. 14) into eq. (IV. 1) gives 


4tr j 





IV. 15 


which has a depolarization component in the plane of incidence in that the 
factor of Jc ^ is not zero; however under the approximation ii = It this result 
readily reduces to eq. (IV. 13) for the case of perfect conductivity. 
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Comparison of the acoustic results given by eq. (IV. 7) with 
eq. (IV. 13) shows that a simulation is possible for the case of back- 
scattering at normal incidence from an imperfectly conducting surface 
under the approximation n = ic . Under the same approximation, it is 
seen that the re-radiation from a perfectly conducting surface can be 
simulated but as before only in- the plane of incidence because of 
depolarization effects. In considering the possibility of simulation out 
of the plane of incidence for the case of a perfect conductor, the same 
conclusion is reached as earlier. However in this case, the quantity 
which has no comparable acoustic term, ( n • !*)]£,, vanishes only upon 
approximation and not as cos© 

OX 

The incident field for the case of polarization parallel to the y-z 
plane is 



A.XIA 


vc 


E.d 




P 






where rr is the same as earlier as are n" , IT , and jT , the unit vectors 
of the local coordinate system. The re-radiated magnetic field which 
contributes to the power flow in the direction of r^ is determined from 





Substituting eq. (IV. 1) into this expression gives the result 

H(P)= f f 'f ^ 


-4TP 


J L 



. IV. 16 
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The factor of r Q in eq. (IV. 14) is neglected here as it was in eq. (IV. 1). 
Making the same severe approximation 


t = i 


P = j 


n = k 


that led to eq. 's (IV. 9) results in 










Using these in eq. 's (IV. 6) and (IV. 5) gives the tangential electric and 
magnetic fields at the surface; these are 



[V. 17a 


T\X H =■ 



[V. 17b 


Substitution of eq. 's (IV. 17) into eq. (IV. 16) gives, again neglecting the 
factor of ~r 

o 


4-rr^, 


X 





><2 




<±S 
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+-(£ r T>) £ 0 -rjl K ' r ds 


IV. 18 

Comparison of eq. (IV. 18) with eq. (IV. 7) leads to the same conclusion 
that was made in connection with the comparison of this equation and 
eq. (IV. 10); i. e. , the similarity exists exactly only in the plane of 
incidence where no depolarization is present and the electromagnetic 
re-radiation process appears to be scalar. The degree of approximation 
in simulation out of the plane of incidence is the same in both cases. 

The approximation made in obtaining the tangential electric 
and magnetic fields at the surface is just as severe in this case of 
parallel polarization as it was earlier. The calculation of the re- 
radiation in the plane of incidence made earlier is not done for this 
case; however, the conclusion made with respect to simulation of normal 
incidence backs cattering is still valid since at normal incidence there 
is no difference in polarization. It remains to determine the magnetic 
field re-radiated from a perfect conductor for the case of parallel polarization. 
Using 


n x E = 0 

nxH=2 (nxIL) 



IV. 19 
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which differs in form with the corresponding result for normal polarization 
given by eq. (IV. 15). A result of these differences is that there is no 
depolarization of the parallel polarized incident wave in the plane of 
incidence (r^ • = 0 ) while, it is recalled, there is for the case of 

normal polarization (however, under the approximation n = It the depolariza- 
tion disappears). These differences in form appear again and must be 
taken into consideration in the simulation of the re-radiation of electro- 
magnetic waves from perfectly conducting surfaces. Comparison of 
eq. 's (IV. 15) and IV. 19) with eq. (IV. 17) shows that simulation of a 
parallel polarized wave requires a perfectly rigid surface ( J?= 1) while 
an inelastic or pressure release surface ( j? = -1) is required for the 
normally polarized case. 

The use of the Kirchhoff approximation has enabled a determination 
to be made of the conditions and the degree of approximation under which 
an acoustic simulation can be made of the re-radiation of electromagnetic 
waves from rough, homogeneous surfaces. Comparison of the analogous 
acoustic and electromagnetic expressions derived under this approximation 
shows in what way similarity exists and the manner in which the vector 
nature of the electromagnetic re-radiation creates the difference existing 
between these two phenomena. For the case of a perfectly conducting 
surface on which both local polarization components are re-radiated equally 
and for the case in which one of the local polarization components was 
ignored (smooth surface approximation) it was found that in the plane of 
incidence the electromagnetic re-radiation is essentially scalar (the 
minor exception to be noted is in connection with eq. (IV. 15)), and an 
acoustic simulation is possible. Out of the plane of incidence for these 
two cases a change in polarization and a difference in form from the 
acoustic occur, but it is still possible to make an approximate simulation. 
From the general investigation of re-radiation in the plane of incidence 
from an arbitrary (yet isotropic) homogeneous surface, it was found that 
the vector re-radiation process creates dependencies on both reflection 
coefficients, and X)/ / which shows that even in the plane of 
incidence, the re-radiation is in general a vector process. However, the 
dependence of the polarized component on the reflection coefficient of 
its orthogonal component is only through the slope dependent terms which 
are also the only terms contributing to depolarization. For backscatter 
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at normal incidence for the general case, there is no dependence upon 
slope terms under the approximation rT= 1 1 with the result that the 
vector nature of the re-radiation disappears and an acoustic simulation 
is possible. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 


A theory ha; 

waves from rough surfaces of the type for which the method of physical 
optics is applicable. The theory was developed for surfaces which are 
generated by stationary, ergodic random processes and which have a 
Gaussian distribution of surface heights. The surfaces are allowed to 
be imperfectly reflecting but are restricted to be homogeneous with 
no variation, random or otherwise, of the material parameters. The 
far-zone power density re-radiated omnidirectionally from the rough 
surface was calculated from the far zone pressure field and the separation 
made between reflected and scattered power densities. The average 
differential scattering cross section, cr Q , was calculated from the 
scattered power density under the condition that the surface was suffi- 
ciently rough to neglect the reflected power. <7 was then determined 
for the cases of Gaussian and exponential dependence of the auto- 
correlation function, and the results of these calculations were seen to 
have strong differences in variation which demonstrates the heavy 
dependence of a on the functional form of the autocorrelation coefficient. 
The expressions for <r 0 are composed of a sum of terms all but one of 
which, the quasi-specular term, have their origin in the first partial 
derivative of the surface under the approximation made. These slope- 
derived terms depend upon the rate of variation of the reflection coefficient 
with angle of incidence, whereas the quasi-specular term depends upon the 
reflection coefficient alone. In and near the specular direction, the quasi- 
specular term was seen to dominate; in other directions, the scattering • 
depends more on the slope-derived terms and consequently on the rate 
of variation of the reflection coefficient with angle of incidence. 

The theoretical values of have been compared with the results 
of experiments that were specially designed to allow a determination of 
the usefulness of the Kirchhoff method. To make the comparison more 
valid, the measured statistical parameters of the surface were used in the 
calculation of the theoretical results. The comparison was made using 
results for <7 which were calculated for Gaussian and exponential 
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autocorrelation functions which bracketed the true autocorrelation co- 
efficient. The evidence obtained by proceeding in this indirect manner 
is sufficiently strong to -justify making the conclusion that the theory is 
a valid one. 

In making the comparison of theory with experiment, the closeness 
of the agreement between the results for the exponential autocorrelation 
function and the experimental data at backscatter is to be noted. If the 
experimental observations had been restricted to backscatter., it would 
be tempting to conclude that the surface is described by an exponential' 
autocorrelation function. However, the poorness of the fit given by the 
use of this function in the specular direction shows that this is not the 
case as do the measurements performed on the rough surfaces. This 
illustrates .the necessity of making omnidirectional measurements. It is 
important to note that to achieve this fit at backscatter it was not necessary 
to use "unreasonable " values for the parameters in the exponential func- 
tion-, i.e. , the parameters used- do describe the surface with some accuracy. 

To determine the similarity of acoustic and electromagnetic re- 
radiation from surfaces, comparison was made between analogous expres- 
sions for the two cases. The electromagnetic expressions were obtained 
by an application of the Kirchhoff method which is a derivation, similar to 
that used to arrive at the acoustic results. It was found that in some 
cases the vector nature of the electromagnetic re-radiation does- not - 
appear and for these cases ah exact acoustic simulation is possible. In 
some other cases where vector effects such, as depolarization and depend- 
ence upon both reflection coefficients, •G and -Q/ ' are present, it 
is still possible, to make an approximate simulation. 
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CHAPTER VI 


RECOMMENDATIONS FOR FURTHER WORK 


In making this study many unsolved rough surface scattering 
problems common to both acoustic and electromagnetics were found. 
The most important of these are the following: 

First and foremost is the development of an approach that will 
yield solutions, even approximate ones, for scattering from surfaces 
of arbitrary roughness. At present, there is no way to treat surfaces 
with irregularities neither small nor large with respect to wavelength. 

An experimental study of the frequency dependence of rough 
surface scattering is badly needed as there has been very little 
verification of the theoretical predictions in hand. A study like this is 
best accomplished acoustically because of the ease with which 
experimental parameters can be varied. 

There- is need for an experimental and theoretical investigation 
of the effect of a layer of different material beneath a rough surface. 
This is another case in which an experimental acoustical study would 
be especially applicable; the experimental results of this work provide 
a starting point for this study. 
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APPENDIX I 

DEFINITION OF CONSTANTS 

The coefficients used in obtaining eq. (II. 14) are defined as 
follows: 

a l - < k lz + k oz> + ^(“» k o 2 - k lz>' 

a 2 “ [ B lx (k oz " k lz> + •f< a « k lx - k ox> - < k lx + k ox>] 

a 3 = [ B ly (k oz - k lz> + X’<“X k ly - V - .ftly + k oy>] 
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APPENDIX II 


CALCULATION OF THE MEANS OF TERMS 
CONTRIBUTING TO THE AVERAGE OF THE RE- RADIATED POWER 


In obtaining the ensemble average of eq. (II. 14), use is first 
made of the ergodicity assumption and, following this, the mean of 
the terms in the integrand is taken. This is done by making use of the 
Karhunen-Loeve theorem CKarhunen, 19473 and the method of character- 
istic functions C Davenport and Root, 1958D . This calculation was first 
made by HoffmanCl955bJ who was one of the first to make a definitive 
study of the statistical difficulties of the scattering problem. 

The Karhunen-Loeve theorem states that a random process 
z = < X 'X) which is continuous in the mean over a finite region D has 

a representation 

Z. = Jl. rvu A II. 1 

where <Z mn Z pq > = 6 mp 6 nq and < Z mn > = 0 if and only if the 
X mn are the eigenvalues and the 0 mn (x,y).the orthonormal eigenfunctions 
of the integral equation 

0()$y) = R & (XjyO 

D 

2 

Here R = r<r is the un-normalized autocorrelation function. It is noted 
that R is a function of both x, y, and x', y 1 and stationarity is not 
assumed at this point. 

This theorem enables the functions of the integrand of eq. (11.14) 
to be written in terms of the $ mn with the properties of the random process 
being given by the expansion coefficients, the Z mn . As an example, 
consider 
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This is written as 


Z*£!. 


AJ 




- i/g__ . 

==■ > 7 w <axp 

3x 




Z ^Z. 




■Z* k-uo (X,y) /\im(a 2 - wir , 


WjV\ 


The average of this quantity is given by 


0 wvi CX>y) Awn Z-Minj^ 

cPx U "~* 


\-l/s 


vMjH 


A II. 2 


To carry the calculation further, the method of characteristic 
functions is introduced. The characteristic function M[jF(x,y)J of the 
function F(x,y) with respect to the random variable Z mn is defined to be 

m[f(x /Y )] = exp 


Therefore , it follows that 
X> ^2.Z-\ 


<*,/)] 

and consequently 


m.n 


< ZiMin e / ~~ l r \y -1 

^ C*»y )J 


Assuming the random process z = X (x,y) is Gaussian with mean zero, 

J - 1/2 

the characteristic function of the function K \ ' <h (x,y) is 

z ran r mn 




exp $\wa (x,y) 


In obtaining this, the relation 
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< Z - /z Y> 2 = 1 

\ mn \ mn// 


which is obtained from above, is used. Therefore 


e' 0 - IT f- ^zT'hMK 1 

m >* iC 4 


<a* 


P f- 7" (x,y) 




ZL, 


The factor of the exponent 



. | j ^ — 


is seen to be 


\A ^ -V- 


<^z- *\ /* /*~ ^P<g ^^vwvi 0<jy)<f> p£ 6*;/ ) 2.wih 2 p£)> 

wtn a<2 




- z 7- foyty pa cxj/j s W p Uj 


vm ) m p, 



which is the variance 


2 

cr 


(x,y)(not a constant) of the process. 


Therefore 

<e^ 2 >=e 


Ki: <rK*,y)/^ 
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and finally 

— x, wi ^ C fyy) e 




Substituting this result into eq. (A II. 2) gives 

' Q 4\ 

ax 






A factor of this is recognized to be the first partial derivative of the 
variance with respect to x ; i. e. 


2jhm« (Vf) g ^ wC> ^y) ^ - a 

IM.iq SX 


So finally 



<9x 



2 

Which, for a stationary process, is zero because cr (x, y) is a constant. 

The evaluation of the means of the terms in the integrand of 

eq. (II. 14) is accomplished by following procedures similar to that used 

above; the calculations are lengthier, however, due to the additional 

complexity of the functions to be averaged. For a stationary, Gaussian 

2 

process with variance cr and normalized autocorrelation coefficient r , 
the results are, using a and p to represent x or y or both 

l<: 2_ ( 

1<Z - — -Z exp 


<^ e 
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<^ze'e '° krif2_2 ' - X, i^cr*- Sr ey p [_ ic?-<r 3 -0 -r)l 
r 32' L 

^ z -( 2 - -2 - N = (~cr a -<35l_ + K£'cr 4 -Sr3r’| 

r L 6>oOp' . aoi.9p-J 

«yp[-i<£cr 3 -(W)] 
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APPENDIX III 


CALCULATION OF THE ORDER OF MAGNITUDE OF 
THE PARTIALLY INTEGRATED TERMS OCCURRING IN THE 
EXPRESSION FOR RE-RADIATED POWER 


The partially integrated terms appearing in eq. (11.21) are of 
two types; these are 



and there are similar terms in which the first integration is performed 
with respect to u . 

The variations of the autocorrelation coefficient and its partial 
derivative in the integrands of the above integrals are given by 


rfay-tX )j rfyy -Jt ) 3 

9 ^ 



9r 



By the assumptions made earlier that the autocorrelation coefficient and 
its partial derivatives decay rapidly relative to the dimensions of the 
aperture, it is clear that these quantities are sensibly zero when the 
variation with respect to y is neglected. Making use of this in 
eq. (A HI. 1) yields the result 
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I 4 [ SXp-f- (ZxM + ^)]j 


-J -J( K-J? 


— ( 6 sw a ^ sivV^Kx^ g 

'Y 


- k£o-z- 


y< 

dxt cjxd 


'/ 


A> Kx~kTy 


A III. 3 


and the integral A III. 2 is zero. 

The result A III. 3 is, neglecting the constants involving the 
material parameters , within a factor of l/K^ of equalling the specular 
power density given by eq. (II. 17). It is because of this factor of l/K^. 
that A III. 3 goes to zero in directions for which = 0 (one of which is 
the specular direction) and the specular power density does not. There- 
fore, the variation of A III. 3 is not the same as the specular power density 
due to the factor 1/K y and also due to the constants involving the 
material parameters, but both vanish in the same way as the surface . 
roughness increases. It is important to note that the result A III. 3 
is similar to the "edge effect" terms found by Beckman 1963 . 

Similar results are obtained for the terms in which the partial 
integrations are first performed with respect to u . 
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APPENDIX IV 

CALCULATION OF THE APERTURE EFFECT 


The aperture effect is evaluated by calculating the integral 
of eq. (III. 4) 

f f ‘grC^g-Cft-) j A 

s r r 3 ' p r~ A IV. l 

for the measured laboratory parameters as a function of the set of angles 
o' ^ 1 ' $ 2 } ’ 

The normalized gain functions, g^,( >3> T ) and # are given by 

the square of the normalized antenna voltage patterns shown in figure 
III. 10. Because of the close fit to the experimental curves obtained with 
the theoretical pattern, the gain functions are, to a very good approximation 




A IV. 2a 



(JCQ 

_ kq 


i ^ — 


A IV. 2b 


where ka = 34 . Substitution of eq. ’s (A IV. 2) into (A IV. 1) and use of 
the coordinate system of figure III. 1 to define the integration gives the 


result 


00 

If 


. kq JSiH A f ? T 


3-r 


Q Q (<q S \ VU 

L kqSt^r J 


— CO 
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The illuminated area S is usually defined, for reasons of convenience, 
on the half power points of the antenna beam and this gives a finite 
illuminated area. It is not necessary to do this here and the limits of 
integration are infinite. The quantities sin^, simE^ , r^, and r f are 
functions of (x,y), the coordinates of the point on the mean plane. 
These functions are 

r T 2 - <x To - x) 2 + <y To - y) 2 + z To 2 

r r 2 “ (x ro - X)Z + ( Vro ' ^ + z ro 2 

2 2 

sin = 1 - cos 

2 2 

sin -\!> r = 1 - cos ^ 


where 


2 2 2 
2. r To +r T " r 

COS = s 

T 2r m r n 


To T 


2 2 2 
o r ^ + r„ - r 

z _ ro r 

cos if = — r 

r 2r r 

ro r 


, 2 2,2 

and r = x + y 

The rectangular coordinates giving the locations of the transmitting and 
receiving antennas are converted from the spherical co-ordinates (r^ 
r ro ' $o ' ^1 ' $ 2 ^ ^ the transformation 

= 0 (principal ray of the transmitting antenna chosen to lie in the 
y-z plane) 

y To ~ r To cos< ^ 1 


z To r To Sin ^ 1 
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X ro = r ro sin ^2 sin ^o 
y ro = r ro sin ^2 cos< ^o 
z ro = r ro cos <» 2 


The integral proved to be too complex to be done analytically and 
it was necessary to resort to a numerical technique of integration which 
was done by machine calculation. The technique is based on the method 
of Gaussian quadratures and is described in the book by Krylov [1962] . 
The results of the machine integration are top lengthy to be included 
here; however, a similar calculation has been performed by the workers 
at the Ohio State University and is available in report form CBarrick , 
1964] . The results are not the same due to a difference in antenna 
patterns; however, the variation as a function of (<f> Q , is 

•similar and representative. 



139 


APPENDIX V 


, DERIVATION OF A VECTOR FORM 
OF THE HELMHOLTZ INTEGRAL 


Stratton [1941] has shown that the electric .field E(P) at- an 
observation point P (see figure IV. 1) due to an electric field IT and 
a magnetic field H distributed over an area S enclosed by a contour 
C is 


Ho) 


t 



s ' 




A V. 1 

where the time variation has been changed from e to e ia)t The 

Green’s function </> for this time variation is 


/ — .'C/t. 2 

ft = ■ A V. 2 

as before and 

V d = ^(<e 

' A V. 3 


approximately for the far field case, k>>^ . In eq. (A V. 3) R^ , is the 
unit vector in the R direction. Stratton has further shown that a surface 
current 


K = -n x H 


A V. 4 
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and a surface charge 

r> = -e x n • E A V. 5 

can be introduced and used; these two quantities are related through the 
continuity equation 

-iwrj — V • K A V. 6 

Combining eq. ‘s (A V. 6) and (A V. 5) and then using eq. A V. 4 gives 

■fi. E. = - V> Cnx'H) 

Aj CO 

Using this and eq. (A V. 3), the last integral on the right hand side of 
eq. (A V. 1) becomes 


ke, 


- J_ nv e 

4i? J ) 4-TPuje, ' ' r 

vS i S - 


xjk.R- . _ — u \ 

V«(nxH)dS 


A V. 7 


Using the identity 

[V. (nxT?)]er° kR _ 

^ R. 


]- 




R. 


in eq. (A V. 7) yields 

i rr. ^ a 


41? JJ 

s 




ke> f fy. Rftxffler 

4-1? cue, J J ^ R. 

I c 


,-^kRI 


dS 


- kj\ f f/ok dS 

41T0U6, 3 ^ 
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Reducing the first integral on the right hand side using Gauss' Law in 
:wo dimensions, this becomes 
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where n ' is the normal unit vector pointing outward from C. 
Since 

n x H • n ' = nx n • H 
and since 

n 1 x n = - <31 
the result is finally 
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A V. 8 

Substituting eq. 's (A V. 8), (A V. 2) and (A V. 3) into eq. (A V. 1) gives 
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A V. 9 


Making S an aperture of the size discussed in Section II. B and manipulating, 
eq. (A V. 9) becomes 



A V. 10 
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